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ABSTRACT 
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rithm for  computing  the  furthest-site  geodesic  Voronoi  diagram  of  k 
point  sites  within  a  simple  n-sided  polygon. 
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1.  Introduction 

A  common  goal  of  much  recent  research  in  computational  geometry  is  to  extend 
algorithms  that  have  been  developed  for  the  Euclidean  metric  to  the  more  compli- 
cated geodesic  metric  inside  a  simple  polygon.  The  geodesic  distance  between  two 
points  in  a  simple  polygon  is  the  length  of  the  shortest  path  connecting  the  points  that 
remains  inside  the  polygon.  For  example,  Toussaint  [T86]  gives  an  algorithm  for  the 
"relative  convex  hull"  of  a  set  of  points  inside  a  simple  polygon;  Aronov  [A87]  gives 
an  algorithm  for  the  nearest-neighbor  geodesic  Voronoi  diagram;  and  Pollack,  Sharir 
and  Rote  [PSR87]  give  an  algorithm  for  the  "geodesic  center"  of  a  simple  polygon. 

A  classic  structure  in  the  Euclidean  metric  is  the  "furthest-site  Voronoi 
diagram."  Given  a  finite  collection  of  point  sites  in  the  plane,  the  furthest-site  Voro- 
noi diagram  partitions  the  plane  into  Voronoi  cells,  one  cell  per  site.  The  site  that 
owns  a  cell  is  the  site  that  is  furthest  from  every  point  in  the  cell.  Using  well-known 
algorithms,  the  Euclidean  furthest-site  Voronoi  diagram  of  k  sites  can  be  computed  in 
time  0(*logjt)  and  space  0{k)  [PS85]. 

The  content  of  this  paper  is  an  efficient  algorithm  for  computing  the  furthest-site 
Voronoi  diagram,  defined  by  the  geodesic  metric  inside  a  simple  polygon.  The  algo- 
rithm uses  0((n+Jk)log(rt+A:))  time  and  0(n+k)  space,  where  n  is  the  number  of 
bounding  edges  of  the  polygon  and  k  is  the  number  of  sites.  The  best  previous  algo- 
rithm for  this  problem  had  running  time  C>(/i-'loglogn)  [AT86],  and  just  computed  (a 
superset  of)  the  vertices  of  the  furthest-site  Voronoi  diagram  of  the  n  corners  of  the 
polygon.  We  remark  that  our  furthest-site  geodesic  Voronoi  diagram  algorithm  is  a 
factor  of  O(logn)  faster  than  the  best  known  nearest-site  geodesic  Voronoi  diagram 
algorithm  [A 87]. 


(')  A  preliminary  version  of  this  paper  appeared  in  Proceedings  of  the  Fourth  ACM  Symposium 
on  Computational  Geometry,  1988. 

(••)  The  work  of  the  first  author  was  supported  by  an  AT&T  Bell  Laboratories  Ph.D.  Scholar- 
ship. Part  of  the  work  was  performed  while  the  first  author  was  at  AT&T  Bell  Laboratories, 
Murray  Hill,  NJ. 


The  problem  of  computing  the  furthest-site  Voronoi  diagram  is  an  extension  of 
the  "furthest  neighbor  problem",  which  is  "Given  a  finite  collection  of  points,  for 
each  point  identify  the  element  in  the  coUeaion  that  is  maximally  distant  from  it." 
Suri  [S87]  shows  how  to  solve  a  special  case  of  the  furthest  neighbor  problem  in  the 
geodesic  metric  inside  a  simple  polygon.  Specifically,  he  gives  an  algorithm  that  for 
each  comer  of  the  polygon  computes  the  corner  that  is  maximally  distant  from  it  in 
the  geodesic  metric.  His  algorithm  runs  in  time  0(n\ogn)  and  space  0(n),  where  n 
is  the  number  of  bounding  edges  of  the  polygon. 

The  geodesic  furthest-site  Voronoi  diagram  generalizes  the  geodesic  furthest- 
neighbor  mapping  of  Suri  [S87]  in  two  ways.  First,  the  Voronoi  diagram  provides  a 
planar  partition  of  the  polygon  together  with  its  interior  into  furthest-site  Voronoi 
cells.  Consequently,  arbitrary  furthest-site  queries  can  be  answered  using  a  planar 
point- location  algorithm.  Second,  the  set  of  sites  is  not  restricted  to  the  comers  of 
the  polygon.  Rather,  the  sites  can  be  arbitrarily  situated  m  the  polygon.  Both  of 
these  generalizations  have  substantial  technical  impact  on  the  algorithm  for  comput- 
ing furthest-site  Voronoi  diagrams. 

There  are  many  analogies  between  the  Euclidean  furthest-site  Voronoi  diagram 
and  the  geodesic  furthest-site  Voronoi  diagram.  In  the  Euclidean  case,  if  a  site  has 
nonempty  Voronoi  cell,  then  it  lies  on  the  convex  hull  of  the  set  of  sites  and  it  does 
not  appear  on  the  line  segment  between  two  sites.  The  counterclockwise  sequence  of 
Voronoi  cells  (at  infinity)  is  the  same  as  the  counterclockwise  sequence  of  sites  on  the 
convex  hull.  In  the  geodesic  case,  we  show  that  a  site  with  nonempty  Voronoi  cell 
lies  on  the  relative  convex  hull  of  the  set  of  sites  and  does  not  appear  on  a  geodesic 
between  two  sites.  The  counterclockwise  order  of  Voronoi  cells  along  the  boundary 
of  the  polygon  is  a  suborder  of  the  counterclockwise  order  of  sites  on  the  relative 
convex  hull.  In  the  Euclidean  case  this  characterization  is  exact:  any  site  on  the  con- 
vex hull  not  between  two  other  sites  has  nonempty  Voronoi  cell.  In  the  geodesic  case 
the  characterization  is  not  exact,  roughly  because  the  polygon  may  not  be  not  large 
enough  for  the  cell  to  appear. 

A  further  analogy  between  the  two  cases  is  the  structure  of  the  Voronoi  diagram 
itself.  In  the  Euclidean  case  the  Voronoi  diagram  forms  a  tree  with  root  at  the 
Euclidean  center  of  the  set  of  sites.  (The  center  of  a  set  of  point  sites  is  the  point 
that  minimizes  the  maximum  distance  to  any  site.)  If  edges  are  directed  towards  the 
root,  then  this  orientation  is  consistent  with  geometric  direction  towards  the  center. 
In  the  geodesic  case  exactly  the  same  properties  hold,  substituting  "geodesic  center" 
for  "Euclidean  center"  and  "geodesic  direction"  for  "direction." 

The  algorithm  for  computing  the  furthest-site  Voronoi  diagram  consists  of  two 
steps.  First,  we  compute  the  restriction  of  the  Voronoi  diagram  to  the  boundary  of 
the  polygon.  Intuitively,  the  boundary  of  the  polygon  in  the  geodesic  case 
corresponds  to  points  "at  infinity"  in  the  Euclidean  case.  Second,  we  extend  the 
diagram  to  the  interior  of  the  polygon.  Because  the  Voronoi  diagram  forms  a  tree 
with  root  at  the  geodesic  center,  the  second  step  is  easy.  It  can  be  performed  by  a 
"reverse  geodesic  sweep"  towards  the  geodesic  center. 
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The  first  step,  the  computation  of  the  Voronoi  diagram  on  the  boundary  of  the 
polygon,  is  much  more  involved.  We  use  a  technique  developed  by  Suri  [S87]  for 
determining  furthest  neighbors.  We  reduce  the  problem  to  three  instances  of  the 
"two-fragment  problem";  an  instance  of  the  two-fragment  problem  consists  of  a  frag- 
ment of  the  boundary  of  the  polygon  and  a  fragment  of  the  relative  convex  hull  of 
the  set  of  sites.  The  relative  convex  hull  fragment  contains  the  furthest  sites  of  all 
points  on  the  polygon  boundary  fragment.  We  solve  an  instance  of  the  two-fragment 
problem  using  divide  and  conquer  in  the  following  manner.  The  polygon  boundary 
fragment  is  split  at  its  midpoint;  this  implies  a  corresponding  split  of  the  convex  hull 
fragment.  Thus  one  instance  of  the  two-fragment  problem  results  in  two  simpler 
instances.  Eventually  instances  become  small  enough  to  be  solved  directly. 

We  refine  Suri's  two-fragment  technique  in  two  ways.  First,  Suri's  algorithm 
always  splits  the  polygon  boundary  fragment  at  a  corner  of  the  polygon.  This  is  suf- 
ficient for  the  furthest-neighbor  problem,  because  furthest-neighbor  information  for 
points  along  a  wall  is  not  of  interest.  For  the  furthest-site  Voronoi  diagram,  simply 
splitting  at  comers  is  insufficient,  since  potentially  many  Voronoi  cells  meet  a  single 
wall  of  the  polygon.  If  necessary,  we  further  split  each  wall  into  subsegments  so  that 
the  shortest  path  tree  from  a  point  in  a  subsegment  to  the  sites  is  combinatorially 
invariant  over  the  entire  subsegment.  The  combinatorial  invariance  of  the  shortest 
path  tree  implies  that  the  Voronoi  partition  of  the  subsegment  can  be  easily  com- 
puted. 

The  second  refinement  of  Suri's  technique  concerns  the  complexity  analysis  of 
the  recursion.  Suri's  original  algorithm  required  a  step  called  "trimming";  trimming  a 
two-fragment  instance  introduces  a  different  subproblem  that  could  be  solved 
directly.  This  operation  is  necessary  in  Suri's  analysis  in  order  to  maintain  the  linear- 
ity of  the  total  size  of  all  subproblems  at  a  particular  level  of  recursion.  We  show 
that  even  without  trimming,  the  total  size  of  all  subproblems  at  a  particular  level  is 
linear.  This  observation  simplifies  the  recursive  structure  of  the  "two-fragment" 
algorithm  so  that  it  actually  matches  the  description  given  above.  The  analysis  has 
been  incorporated  into  Suri's  furthest-neighbor  algorithm. 

The  best  lower  bound  that  we  know  for  computing  the  furthest-site  geodesic 
Voronoi  diagram  is  Oin+k\ogk).  It  follows  from  known  lower  bounds  for  diameter 
computation  in  the  Euclidean  case.  Conceivably,  the  current  algorithm  could  be 
improved  to  match  this  lower  bound. 

2.  The  Furthest-Point  Geodesic  Voronoi  Diagram 

This  Section  contains  the  definition  of  the  furthest-point  geodesic  Voronoi 
diagram  and  some  of  its  basic  properties.  We  begin  by  discussing  geodesies  in  Sec- 
tion 2.2.  Section  2.3  contains  a  fairly  extensive  treatment  of  relative  convex  hulls. 
The  notion  of  a  "far  side"  of  a  relative  convex  hull  is  developed  in  Section  2.4;  this  is 
a  technical  idea  used  to  prove  the  Ordering  Lemma  in  Section  2.7.  Sites  are  assumed 
to  be  in  "general  position";  this  assumption  and  its  consequences  are  discussed  in  Sec- 
tion 2.5.    In  Section  2.6  we  actually  give  the  definition  of  the  Voronoi  diagram. 


Section  2.7  contains  the  Ordering  Lemma,  which  states  that  the  order  of  Voronoi 
cells  around  the  boundary  of  the  containing  polygon  is  the  same  as  the  order  of  sites 
around  the  relative  convex  hull  of  the  set  of  sites.  In  Section  2.8  we  define  a  refined 
form  of  the  Voronoi  diagram  and  use  it  to  show  a  linear  bound  on  the  descriptive 
complexity  of  the  (unrefined)  Voronoi  diagram.  Finally,  in  Section  2.9  we  show  that 
the  Voronoi  diagram  forms  a  tree  directed  towards  the  geodesic  center  of  the  set  of 
sites.  The  algorithm  for  computing  the  diagram  appears  in  Section  3. 

2.1.  Preliminaries 

The  universe  U  is  a.  compact  region  in  the  plane  whose  boundary  dU  is  a.  simple 
n-sided  polygon.  S,  the  set  of  sites,  is  a  collection  of  k  points  of  U.  A  vertex  of  dU 
is  called  a  comer  and  a  segment  of  3  C/  is  a  wall.  A  comer  is  reflex  if  the  measure  of 
its  interior  angle  is  more  than  it  and  convex  if  it  is  less  than  ir.  If  jc  and  y  are  distinct 
points  of  dU,  then  boundary  fragment  dU[x,y]  is  the  portion  of  dU  counterclockwise 
from  x  to  y  inclusive.  The  symbol  d  denotes  the  boundary  of  a  set  relative  to  the 
whole  plane,  rather  than  to  any  proper  subset  of  the  plane.  The  terms  "relative 
boundary"  and  "relative  interior"  used  without  any  other  qualification  mean  "relative 
to  U." 

A  polygonal  path  is  a  simple  path  comprised  of  a  sequence  of  line  segments.  If 
p  is  a  polygonal  path,  then  the  size  of  p,  |p  |,  is  the  number  of  maximal  segments 
contained  in  p  and  not  containing  a  comer  of  dUin  their  interiors,  and  the  length  of  p 
is  the  sum  of  their  (Euclidean)  lengths.  A  polygonal  region  is  a  compact  set  whose 
boundary  is  the  union  of  a  finite  number  of  line  segments.  We  allow  points  as 
(degenerate)  line  segments,  so  for  example  a  finite  set  of  points  is  also  a  polygonal 
region. 

2.2.  Geodesies 

For  points  x,y^U,  the  geodesic  path  g(x,y)  is  the  shortest  path  in  C/ connecting  x 
and  y.  Such  a  shortest  path  is  called  simply  a  geodesic.  This  path  is  unique.  In  fact, 
it  is  a  polygonal  path  with  interior  vertices  only  at  reflex  comers  of  dU  [LP84].  We 
often  consider  g(x,y)  directed  from  x  to  y.  A  link  of  g(x,y)  is  a  maximal  segment  of 
gix,y)  not  containing  any  comers  in  its  interior;  clearly  an  endpoint  of  a  link  is  either 
X,  y,  or  a  comer.  The  first  endpoint  of  the  last  link  of  ^(x,y)  is  the  anchor  of  y  with 
respect  to  x;  it  is  either  a  reflex  comer  of  flf/  or  jc  itself. 

We  make  heavy  use  of  the  fact  that  the  intersection  of  any  two  geodesies  is  con- 
nected (and  is  itself  a  geodesic).  This  fact  follows  immediately  from  uniqueness  of 
geodesies.   Two  geodesies  overlap  if  they  intersect  in  more  than  a  single  point. 

The  geodesic  distance  d(x,y)  between  points  x  and  y  is  the  length  of  gix,y). 
The  geodesic  distance  is  a  metric;  in  particular,  it  is  continuous  as  a  function  of  both  x 
and  y  and  satisfies  the  triangle  inequality  dix,z)^d(x,y)+d(y,z).  Furthermore,  by 
uniqueness  of  geodesies,  dix,y)  +  d(y,z)  =  d(x,z)  if  and  only  if  y  lies  on  g{x,z).  We 
often  write  du  for  the  function  defined  by  du{x)  =  d(u,x). 
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Lemma  2.2.1:  [PSR87,  Lemma  1]  For  any  u,v,w^U,  du  is  a  convex  function  on 
giv,w)  with  unique  local  minimum  (possibly  at  v  or  w).  In  particular,  for  any 
z^giv,w),  z¥'v,w,  duiz)<mix{du(v),du(w)}. 

The  geodesic  direction  B(x,y)  from  x  to  y#jc  is  the  direction  (i.e.,  unit  vector) 
from  X  towards  the  anchor  of  x  with  respect  to  y.  For  fixed  y,  Bix,y)  =  -Vdy,  where 
Vdy  is  the  gradient  of  dy  with  respect  to  x,  evaluated  at  x  [A87,  3.12].  At  any  point  x 
not  an  anchor  with  respect  to  y,  dy  is  differentiable  as  a  function  of  x  and  9(x,y)  is 
continuous  as  a  function  of  both  x  and  y.  For  H  a  closed  subset  of  U,  let  Q{x,H)  be 
{Q(x,h):h^H,  hi^x). 

The  geodesic  angle  Lxyz  is  the  angle  counterclockwise  from  B{y,x)  to  B{y,z).  The 
measure  of  Lxyz  will  be  written  mLxyz.  The  angle  between  B{y,x)  and  6(y,z)  is  the 
smaller  of  the  two  angles  Lxyz  and  Lzyx. 

Lemma  2.2.2:  [PSR87,  Corollary  2]  lix,y,z^U,  y^x,z,  and  the  angle  between  9(y,j:) 
and  8(y,z)  is  at  least  11/2,  then  d(x,z)  >  max{d(x,y),d(y,z)}. 

The  shortest  path  tree  from  s,  T{s),  is  the  union  of  the  sets  of  links  of  gi.s,y) 
taken  over  corners  y  of  dU.  It  has  n  — 1  or  n  links,  depending  on  whether  or  not  s 
itself  is  a  comer  of  U  [GHLST]. 

Let  Pais)  be  the  set  of  points  in  U  that  have  anchor  a, with  respect  to  s.  The 
shortest  path  partition  of  U  from  s  is  the  collection  {Pais):  Pais)¥'0}.  It  is  a  planar 
polygonal  subdivision  of  U;  it  can  be  computed  and  in  fact  triangulated  in  linear  time 
given  a  triangulation  of  U  [GHLST].  We  can  describe  the  bounding  edges  of  the 
shortest  path  partition  as  follows.  Suppose  Pais)  is  not  empty,  where  a #5.  Let  ab 
be  the  first  link  of  the  geodesic  gia,s)  (clearly  ab  is  the  second  link  of  all  geodesies 
gix,s)  for  x(.Pais)).  Since  Pais)  is  not  empty  and  a  is  a  reflex  comer  of  dU,  we  can 
extend  link  ab  past  a  into  U.  First  suppose  that  neither  wall  of  dU  incident  to  a  over- 
laps this  extension.  Let  y  be  the  first  point  past  a  of  the  extension  so  that  y^dU. 
Then  segment  ay  is  the  shortest  path  partition  edge  (from  s  with  anchor  a),  denoted 
Pais).  Now  suppose  some  wall  of  dU  overlaps  the  extension  of  segment  ab;  then  we 
simply  define  pais)  to  be  this  wall.  It  can  be  checked  that  the  boundary  of  a  cell  of 
the  shortest  path  partition  consists  of  an  altemating  sequence  of  shortest  path  parti- 
tion edges  and  sections  of  dU. 

A  set  AQU  is  relatively  convex  (with  respect  to  U)  if  gix,y)CA  whenever  x,y  €A; 
the  relative  convex  hull  of  set  F,  denoted  RiF),  is  the  smallest  relatively  convex  set 
containing  F  (that  is,  the  intersection  of  all  relatively  convex  sets  containing  F)  [T86]. 
Relatively  convex  sets  are  discussed  in  detail  in  the  next  section.  A  set  is  degenerate 
if  it  is  contained  in  a  single  geodesic. 

If  X  and  y  are  distinct  points  of  U,  then  a  shadow  of  gix,y)  is  a  point  y^dU  so 
that  segment  yy  extends  the  last  link  of  g  ix,y)  while  staying  in  U.  Clearly  y  is  a  sha- 
dow of  ^(;c,y)  only  if  y  €317.  Shadows  are  not  unique,  indeed  it  is  possible  that  every 
point  on  a  subsegment  of  a  wall  is  a  shadow  of  gix,y).  Similarly,  a  foreshadow  of 
gix,y)  is  a  point  x'^dU  lying  on  a  segment  contained  in  U  extending  the  first  link  of 
giXyy)  backwards. 
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A  boundary  geodesic  is  a  geodesic  connecting  two  distinct  points  of  dU.  Let  x* 
and  y  be  the  foreshadow  closest  to  x  and  shadow  closest  to  y  of  g(x,y),  respectively. 
Then  g(x,y)  denotes  the  boundary  geodesic  gix'.y).  Geodesic  g(x,y)  splits  U  into 
two  simply  connected  polygonal  regions  U[x,y]  and  U\y,x]  with  disjoint  interiors; 
d(U[x,y])  is  dU[x\y]yjg(y,x')  and  d(U\y,x])  is  g(x' ,y)\JdU\y,x'].  (Note  that 
diU[x,y])  is  distinct  from  dU[x,y].)  Intuitively,  U[x,y]  contains  points  lying  on  or  to 
the  right  of  gix',y),  while  points  on  the  geodesic  or  to  the  left  of  it  constitute 
U\y,x].  Notice  that  ^(;c,>')  is  exaaly  the  common  boundary  of  U[x,y]  and  U[y,x]; 
hence  any  geodesic  from  a  point  in  U[x,y]  to  a  point  in  U[y,x]  must  intersect  gix,y). 
Both  U[x,y]  and  U\y,x]  are  relatively  convex,  since  a  geodesic  connecting  two  points 
of,  say,  U[x,y]  must  have  connected  intersection  with  gix,y). 
Lemma  2.2.3:  Suppose  u,v^dU,  u^v,  wig(u,v),  x^U[u,v],  and  wt^x. 

(1)  If  xig(u,v),  then  any  shadow  3c  of  g(H',ar)  lies  in  dU[u,v]. 

(2)  If  Ar^^(M,v),  then  some  shadow  x  of  ^(>v,jr)  lies  in  dU[u,v]. 

Proof:  (1)  Suppose  xigiu,v).  Then  g{w,x)  cannot  intersect  g(u,v)  again  after  x,  so 
x^dU[u,v]. 

(2)  Suppose  xigiu.v);  without  loss  of  generality  assume  w,  x,  v  are  in  that  order 
along  g(u,v).  We  can  choose  x=v  unless  ^(u,v)  bends  at  or  after  x.  If  ^(m,v)  bends 
right  at  some  point  c  at  or  after  x,  then  since  U[u,v]  lies  locally  to  the  right  of 
^(m,v),  c  must  be  a  reflex  comer  of  dU[u,v],  and  we  can  choose  x=c.  If  g(u,v) 
bends  left  at  some  point  c  then  the  straight-line  continuation  of  g(w,c)  at  c  enters  the 
interior  of  U[u,v]  and  thus  will  not  intersect  giu,v)  again.  Hence  we  can  choose  I  to 
be  any  shadow  of  giw,c)  distinct  from  c.   □ 

Lemma  2.2.4:  Suppose  u,v€3C/,  u^v,  w^g{u,v),  x^U[u,v],  w^x,  and  w  the  closest 
shadow  of  ^(j:,w). 

(1)  If  widU and  g(x,w)  does  not  overlap  giu,v),  then  w^dU[v,u]. 

(2)  If  V  is  the  closest  shadow  of  g(u,w)  and  giu,w)  is  not  an  initial  portion  of  giu,x) 
(i.e.  wigiu,x)),  then  w€dU[v,u]. 

Proof:  (1)  If  gix.w)  does  not  overlap  g(u,v),  then  9(w,x)  cannot  be  8(w,u)  or 
9(w,v).  Since  widU,  vv#w  and  Q(w,w)=  —Q(w,x)  must  enter  the  interior  of  U[v,u] 
at  w.  As  g{x,w)  cannot  intersect  g{u,v)  again,  w^dU[v,u]. 

(2)  The  statement  is  trivial  if  w^dU.  If  not,  then  w¥=w.  We  might  have 
Q(w,x)  =  Q{w,u),  in  which  case  vv  =  v.  We  cannot  have  9(w,x)=  -9(w,«),  else 
w^g(u,x).  Otherwise  9(w,vv)  must  enter  the  interior  of  U[v,u]  and  as  before 
w^dU[v,u].  n 

Boundary  geodesic  g{x,y)  separates  points  a  and  b  if  a^U[x,y]  and  b^U[y,x], 
or  vice  versa.  Similarly,  boundary  geodesic  gix,y)  separates  sets  A  and  B  if 
ACU[x,y]  and  5Cf/Ijy,x],  or  vice  versa.  Note  that  separation  does  not  imply  dis- 
jointness,  indeed  if  (degenerate)  set  A  is  contained  in  boundary  geodesic  g(x,y),  then 
g(x,y)  separates  A  from  itself. 

For  a  compact  set  FCU  and  z^U,  let  ra^(z,F)  =  sup  d^ix).  The  center  of  F  is 
the  point  z  that  minimizes  rad{z,F).   Pollack,  Sharir,  and  Rote[PSR87]  show  that  the 
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center  of  the  set  of  vertices  of  U  is  unique.  In  fact,  with  minor  modifications  their 
proof  shows  that  the  center  of  any  compact  set  FCt/  is  unique. 

For  points  x,y,z€U  we  define  the  geodesic  triangle  Ajcyz  as  follows  [PSR].  If 
{x,y,z}  is  degenerate,  then  Ajcyz  is  the  smallest  geodesic  containing  x,  y,  and  z.  Oth- 
erwise we  can  choose  points  x' ,  y' ,  and  z'  so  that  x'  is  the  point  at  which  g  (x,y)  and 
g(x,z)  diverge,  and  similarly  for  y'  and  z'.  Then  the  circuit  gix',y'),  g(y',z'), 
g(z',x')  is  a  simple  polygon  Ax'y'z'.  We  define  Axyz  to  be  the  union  of  gix,x'), 
^(y.y').  g(z,z')f  and  Ax'y'z'  together  with  its  interior.  All  of  the  interior  angles  of 
Ax'y'z'  are  reflex  except  at  x' ,  y' ,  z'  [PSR].  The  geodesic  triangle  is  a  special  case  of 
the  relative  convex  hull  of  a  finite  set  of  points  (discussed  in  Section  2.3).  The  proof 
of  the  following  Lemma  is  in  the  same  spirit  as  the  proofs  of  Lemmas  2.2.3  and 
2.2.4. 

Lemma  2.2.5  (Triangle  Lemma):  Suppose  {x,y,z}  is  not  degenerate.  Let  y  and  z  be 
shadows  of  g(x,y)  and  gix,z)  respectively.  Assume  ziU[y,x].  If  w^Aryz  and 
uig(x,y)ng{x,z)  then  there  is  a  shadow  u  of  gix,u)  so  that  u^dU]y,'z]  and  g(x,u) 
intersects  giy.z).  If  u  is  in  the  interior  of  Axyz,  then  for  any  shadow  u  of  g(x,u), 
u€dU\y,z]  and  g(x,ii)  intersects  giy,z).  If  u^U[y,x]nU[x,z]nU\y,z],  then  g{x,u) 
interseas  g(y,z)  and  there  is  a  shadow  u  of  g(x,u)  in  dU\y,'z]. 

2.3.  Relatively  Convex  Sets 

This  section  develops  properties  of  relatively  convex  sets.  The  main  result  is 
Lemma  2.3.4.  It  states  that  the  "extreme"  elements  of  a  set  F  can  be  ordered  so  that 
the  relative  convex  hull  of  F  is  the  intersection  of  all  "cones"  defined  by  consecutive 
triples  of  extreme  elements.  An  immediate  consequence  of  Lemma  2.3.4  is  a  decom- 
position of  a  relatively  convex  set  into  a  collection  of  simple  polygons  and  connecting 
geodesies.  Also,  the  order  of  extreme  elements  extends  to  a  natural  notion  of  a 
traversal  of  the  boundary  of  a  relatively  convex  set.  This  ordering  has  a  number  of 
useful  consequences,  given  in  Lemmas  2.3.6  through  2.3.9. 

Lemma  2.3.1:  Any  relatively  convex  set  R  is  simply  connected. 
Proof:  We  show  that  the  region  enclosed  by  a  simple  cycle  y  in  R  lies  completely  in 
R.   Suppose  x^7,  y  is  in  the  interior  of  7,  and  y  is  a  shadow  of  g(x,y).   Since  yy  con- 
nects a  point  inside  7  to  a  point  on  or  outside  7,  there  is  a  point  w  in  the  intersection 
of  yy  and  7.   Since /?  is  relatively  convex,  g(x,w)CR,  so  y^R.   □ 

Let  F  be  a  nonempty  polygonal  region  contained  in  U.  For  x^U,  we  define 
spaniXyF)  to  be  the  smallest  set  of  directions  containing  Q(.x,F)  so  that  whenever 
a,a'€9(x,F),  a^  —a',  and  the  angle  between  a  and  a'  is  contained  in  U  near  x, 
then  every  direction  in  this  angle  is  in  span(x,F).  If  xdU  is  not  a  reflex  corner  of 
dU,  then  it  is  easy  to  see  that  either  span{x,F)  has  a  single  component  or  span(x,F) 
consists  of  two  opposite  directions.  If  j:  is  a  reflex  corner  of  dU,  then  span  (x,F)  may 
have  two  components  that  are  not  opposite  dirertions;  in  fact  one  or  both  components 
may  have  positive  measure.  Since  F  is  a  polygonal  region,  it  is  easy  to  check  that  the 
endpoint  of  any  component  of  span(x,F)  is  9{x,y)  for  some  y€F.  We  show  below 
(Lemma  2.3.3)  that  span (x,F)  =  Qix,R  (F)). 


If  span{x,F)  consists  of  a  single  component  of  measure  less  than  2it,  then  x  is 
an  exterior  point  of  F.  Notice  span  {x,F)  can  have  measure  exceeding  it  but  less  than 
2-11  only  if  X  is  a  reflex  comer  of  dU.  If  span(x,F)  consists  of  a  single  component  of 
measure  less  than  ir  and  x^F  then  ;c  is  an  extreme  point  of  F.  If  span{x,F)  consists 
of  two  connected  components,  then  x  is  a  thin  point  of  F.  In  Figure  2.3.1,  for 
F={1,2,3,4,5},  point  5  and  all  points  of  dU  except  a  are  exterior  (but  none  are 
extreme),  points  1,  2,  and  3  are  extreme  points,  and  every  point  on  g{\,a)  except  1 
is  a  thin  point.  Set  F  is  extreme  if  every  element  of  it  is  extreme.  In  Figure  2.3.1, 
{1,2,  3}  is  extreme. 


Figure  2.3.1.  Relative  convex  hull 

For  the  remainder  of  this  section,  FCC/  is  a  finite  set  of  points  containing  at 
least  two  elements. 

Let  a:  be  an  exterior  point  of  F  or  a  point  of  dU.  If  x^dU  define  oppfix)  to  be 
X,  otherwise  define  opppix)  to  be  the  first  point  of  dU  intersected  by  the  ray  with 
endpoint  x  directed  opposite  the  bisector  of  span{x,F).  The  clockwise  extreme  point 
of  F  from  x,  denoted  r(x),  is  the  point  f^F  so  that  there  is  no  geodesic  extending 
g{x,f)  to  a  point /'€F,  and  among  all  such  points  the  closest  shadow  of  g(x,f)  is  as 
clockwise  as  possible  in  dU—{oppF(x)}.  Similarly  define  the  counterclockwise  extreme 
point  of  F  from  x,  denoted  l{x).  For  example,  in  Figure  2.3.1,  /(2)  =  1,  r(2)  =  3, 
/(a)  =  1,  r(a)  =  2.  There  are  two  subtleties  to  these  definitions.  First,  r  and  / 
depend  upon  the  set  F,  but  except  in  the  proof  of  Lemma  2.4.3,  we  leave  this  depen- 
dence implicit.  Second,  there  are  two  distinct "  (though  overlapping)  cases  in  the 
definition:  either  x  is  exterior  (and  span{x,F)  consists  of  a  single  component)  or 
xidU  (and  even  if  span(x,F)  consists  of  two  components,  there  is  still  a  natural  way 
to  define  r(x)  and  l(x)).  Notice  there  is  no  natural  definition  of  r(x)  and  l(x)  if 
span(x,F)  has  measure  2it  or  if  span(x,F)  has  two  components  and  xidU.  The  fol- 
lowing Lemma  describes  properties  of  functions  r  and  /;  its  proof  is  quite  technical 
and  detailed. 

Lemma  2.3.2:  Let  x  be  an  exterior  point  of  F  or  a  point  of  dU,  r  =  r{x),  l  =  l(x),  r  the 
closest  shadow  of  g(x,r),  and  /  the  closest  shadow  of  g(x,l). 

1.     If    x     is     exterior     then     span(x,F)=Lrxl.      If    xidU    or     mLrxli^Tt     then 
FCC/[r,x]nC/[x,/];  if  mLrxl  =  Tr  then  FCU[r,l]. 
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2.  If  r#/,  then  rr  and  //  are  disjoint;  if  also  x^dUthen  x,  r,  /  are  in  that  counterclock- 
wise order  around  dU. 

3.  Both  r  and  /  are  extreme  points  of  F. 

4.  If  X  is  an  extreme  point  of  F,  then  r(l(x))=l(r(x))=x. 

5.  Let  x€dUoT  mLrxK-a.  Then  /=r(x)  if  and  only  if /€F,  FQU\f,x],  and  g{x,f) 
cannot  be  extended  to  g{x,f')  for  any  other /'€F. 

Proof:  (1)  We  must  always  have  span{x,F)QLrxl  by  definition  of  span(x,F)  and  of 
points  r  and  /;  certainly  Q(x,r),Q(x,l)€span(x,F).  If  a:  is  exterior  then  span{x,F)  is 
connected,  so  we  must  have  Lrxl=span{x,F).  For  the  second  statement,  first  sup- 
pose xidU  or  mLrxli'it.  Let  x*  be  the  closest  foreshadow  of  g{x,r).  We  claim 
l^dU\r,x*\.  this  follows  immediately  from  the  definition  of  r  and  /  if  mLrxK-n, 
x^dU,  or  mLrxl>'n  (since  in  the  last  case  necessarily  x^dU).  For  any  /€F,  either 
f^gix.r),  f^g(x,l),  or  the  closest  shadow  of  g(x,f)  lies  in  dU\r,i]CdU\r,x*],  so 
f^U[r,x],  and  FCU[r,x].  Similarly  FCU[x,l].  The  argument  that  mLrxl='Tr  implies 
FCU[r,l]  is  similar. 

(2)  Suppose  /#r.  If  //and  rr  met,  either  /  would  lie  on  g(x,r)  or  vice  versa,  contrary 
to  the  definition  of  r,  /.  The  ordering  of  x,  7,  and  /  follows  immediately  by  defini- 
tion. 

(3)  We  assume  FCU[r,x],  otherwise  a  similar  argument  works  using  U[r,l].  To 
show  r  extreme,  we  need  to  show  span(r,F)  is  connected  and  has  measure  less  than 
IT.  Whether  or  not  r^dU,  Q{r,U[r,x])  is  connected.  As  FCU[r,x], 
spanir,F)CQ(r,U[r,x]).  Hence  it  suffices  to  show  that  for  f^F,  6(r,/)  lies  in  the 
angle  from  6(r,x)  clockwise  to  but  not  including  —Q(r,x).  Now  9(r,/)  cannot  be 
-9(r,x)  else  g(x,f)  would  extend  g(,x,r).  Also  6(r,/)  can  be  clockwise  of  —Q{r,x) 
only  in  the  case  that  r  is  a  reflex  comer  of  dU  lying  to  the  right  of  g(x,r);  again  this 
is  impossible  because  g{x,f)  would  extend  g(x,r). 

(4)  By  (1),  FC,U[r,x].  For  any  /€F  not  appearing  on  gix.r),  f  is  in  the  relative  inte- 
rior of  U[r,x].  Thus,  by  Lemma  2.2.3(1),  the  closest  shadow  of  g(,r,f)  lies  on 
3C/[r,jr*],  where  x*  is  the  closest  foreshadow  of  g(x,r).  This  implies  x  =  l(r)=l(r(x)) 
by  definition  of  counterclockwise  extreme  point. 

(5)  Clearly  there  can  be  at  most  one  point  /in  F  satisfying  "FCU[f,x]  and  g(x,f) 
cannot  be  extended  to  g(x,f')  for  any  f'^F".  Since  f=r(x)  is  one  such  point,  the 
claim  follows,  n 

Lemma  2.3.3:  For  x^U,  span{x,F)  =  Qix,R(F)). 

Proof:  Suppose  y,z^F  are  such  that  0<mLyxz<T:  and  Lyxz  is  (locally  around  x)  con- 
tained in  U.  By  examining  Aoryz  and  using  g(y,z)CR(F)  we  see  that 
LyxzCQix,R{F)).  Hence  spanix,F)CQ(x,RiF)). 

If  span(x,F)  has  measure  lit,  then  it  must  be  that  span(x,F)  =  Q(x,R(F)).  Sup- 
pose x  is  an  exterior  point  of  F;  let  r=r(x)  and  l=l(x).  If  mLrxl^^TT  then 
RiF)CU[r,x]r\U[x,l]  since  FCU[r,x]nU[x,l]  by  Lemma  2.3.2(1),  and  U[r,x]  and 
U[x,l]  are  relatively  convex.  Hence  Q{x,R(F))  C  Q{x,U[r,x]nU[x,l])  =  Lrxl  = 
span(x,F),  where  the  first  equality  follows  by  definition  and  the  second  by  Lemma 
2.3.2(1).    If  mLrxl  =  tt  then  R(F)CU[r,l],  so  9(x,R(F))   C   Q(x,U[r,l])   =  Lrxl  = 
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span(x,F). 

Finally,  if  x  is  a  thin  point  of  F,  we  can  find  a  segment  through  x  splitting  U  into 
two  simple  polygons,  each  containing  one  of  the  components  of  span(x,F).  Let  Fi 
and  F2  be  the  intersections  of  these  two  polygons  with  F,  and  let  ^i,  H2  be  their 
respective  convex  hulls  (relative  to  U).  It  is  easy  to  verify  that  x  is  an  exterior  point 
of  both  Fi  and  F2  and  thus,  by  the  first  part  of  the  proof, 

span  (x,Fi)  U  span  {x,Fi)  =  %{x,H\)  U  e(x,ff  2)  Q  Q(x.R  (F)) 

The  inclusion  is  proper  only  if  there  are  two  directions,  Qi^  spanix,Fi)  and 
82  €  span{x,F2)  the  angle  between  which  has  measure  less  than  ir  and  is  contained  in 
U  locally  near  x.  But  then  this  angle  would  be  included  in  span(x,F),  contradicting 
the  assumption  that  a:  is  a  thin  point,  i.e.,  span(x,F)  is  disconnected.   Q 

It  is  immediate  that  any  two  sets  with  the  same  relative  convex  hull  have  the 
same  extreme,  exterior,  and  thin  points. 

The  next  Lemma  is  the  main  result  of  this  section.  To  state  it  we  need  the  notion 
of  a  geodesic  cone.  Suppose  mLxyz<'n  and  ziU[x,y].  Let  s  be  the  segment  from  y 
to  the  closest  foreshadow  y*  of  g(y,x),  open  at  y  and  closed  at  y*.  The  geodesic  cone 
U[x,y,z]  is  U[x,y]r\U[y,z]-s.  See  Figure  2.3.2.  Notice  s  intersects  U[x,y]nU\y,z] 
only  if  s  is  also  a  foreshadow  segment  of  giy,z).  It  is  easily  checked  that  U[x,y,z]  is 
relatively  convex  and  that  the  boundary  of  U[x,y,z]  consists  of  g(z,y),  g(y,x),  and 
dUlx,!],  where  x  and  z  are  the  closest  shadows  of  g{y,x)  and  g{y,z),  respectively. 
Furthermore,  if /is  an  extreme  point  of  F,  then  FCU[r(f),f,l(f)]. 


Figure  2.3.2.  C/[2,l,3]  (solid)  and  U[3,2,l]  (dashed) 


Lemma  2.3.4:  H=R(F)  is  a  simply  connected  polygonal  region.    The  extreme  ele- 
ments of  F  can  be  labelled  /i,  h  =  r(Ji),    ...,  fQ=fm  =  r(fm-i),  fi=fm  +  l  =  rifm)  so 

m  m 

that  H=  nu\fi^^,fi,fi-i]  and  dH=Ug(fi,fi  +  i). 
1=1  1=1 

Proof:  F  must  contain  some  extreme  element,  since  r(u)  is  extreme  for  any  u^dU. 

Let  /i  be  an  extreme  element  of  F,  and  consider  the  sequence  f2  =  r(J\),  /3  =  r(/2), 

Since  F  is  finite  and  Ur  is  the  identity  (Lemma  2.3.2(4)),  it  must  be  that 

fm  +  i=^fl  for  some  m  with  /i, ...,/«  distinct.    If  F  is  degenerate,  then  m=2  and  the 

Lemma  follows  easily.    We  assume  F  is  nondegenerate,  hence  m>2  and  points  in 
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each  triple /i-i,/i,/i+i  are  distinct. 


fli 


Let  /=  n  U\fi+i,'fi,fi-i]-   Our  goal  is  to  show  81=  Ug(fi,fi+i).   Let  gifiJi+i) 

i  =  1  1  =  1 

have  closest  foreshadow  j,-  and  closest  shadow  r,-.  Then  3(C/|/i+i,j^-,/i-i])  = 
g(fi,ti)^gisi-i.fi)UdU[ti.Si-u.  By  Lemma  2.3.2(1)  we  have  for  ;  =  l,...,m, 
fj-i^U[fj  +  i,fj]=U[tj,Sj],  so  using  Lemmas  2.2.3  and  2.2.4,  we  have  Sj-\,Sj,tj-i,tj 
in     that    counterclockwise     order     on     dU    (possibly     Sj=tj-i).      Consequently, 

IR  fft  Fit 

U^dU[sj,tj]=dU.    Hence   r\^dU[tj,Sj-i]  C   n^dU[tj.Sj]  Q  {su...,s„,ti,...,t„}.    We 

m 

now   claim    that    dU[ti,Si-i]r\I   C     U  gisj.ti).     Suppose    x^dU[ti,Si-i]r[I;    either 

7  =  1 
m  m 

x^  n  3C/[f.,J.-i],  in  which  case  the  claim  is  immediate,  or  xi  D  dU[tj,Sj-i].   In  the 
;=i  7=1 

second  case,  xidU[tj,Sj-i]  for  some  j,  but  x^U{fj+i,fj,fj-i]  since  x^I,  and  since 
x^dU,  it  must  be  that  x^d(U\fj+i,fj,fj-i]),  so  x€g(fj,tj)\Jg(sj-i,fj),  establishing 

m 

the  claim.   We  now  show  37=  U5(/i,/l+i).  By  elementary  topology, 

m 

di  =  Uind(u\fi+ufi.fi-i]) 
1=1 

m  m 

=  U(g(fi.ti)ug(fi,si-i)nr)  u  Uoc;[f.-.j,-i]n/) 
1=1  j=i 

m 

=  U^(5,-,/.)n/. 
»=i 

Since  segments  stft  and  ri_i/i  are  in  f/[/i+i,/i,.^-i]  only  if  Si=ti-i=fi€dU,  we  have 

3/=.U.?C/i'/;+i)- 
«=i 

For  each  /,  we  have  HCU[fi+i,fi,fi-i]  since  C/[/)  +  i,//,/i-i]  is  relatively  convex 

and  contains  F.   Hence  HQI.   Also  dICH  since  each  geodesic  is  contained  in  H;  since 

H  is  simply  connected,  H=[.   Every  extreme  element  /of  F  must  be  in  H  and  in  fact 

on  dH,  since  /is  on  the  boundary  of  'J[r(f),f,l(f)].    Since  no  extreme  element  can 

lie  on  a  geodesic  between  two  othei  extreme  elements,  each  extreme  element  is  equal 

to  /•  for  some  i.  o 

Since  /?(F)  is  a  connected  polygonal  region,  it  can  be  decomposed  into  plateaus 
and  bridges:  a  plateau  is  a  maxhnal  closed  two-dimensional  region  whose  boundary  is 
a  simple  polygon  and  a  bridge  Is  a  maximal  polygonal  path  not  containing  any  pla- 
teau points  in  its  interior.  We  require  that  both  endpoints  of  a  bridge  be  in  a  plateau, 
thus  a  single  point  at  the  end  of  a  bridge  may  be  a  plateau.  In  Figure  2.3.2,  point  1 
and  geodesic  triangle  A23a  are  plateaus;  segment  la  is  a  bridge.  Each  plateau  is  rela- 
tively convex  and  each  bridge  is  a  geodesic.  Since  R{F)  is  simply  connected,  this 
decomposition  forms  a  tree,  with  plateaus  taken  as  nodes  and  bridges  as  edges.  An 
extreme  point  x  of  R{F)  is  a  convex  vertex  of  a  plateau  (if  span(x,F)  has  positive 
measure)  or  is  a  plateau  itself  (if  span{x,F)  consists  of  a  single  direction).    A  thin 
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point  X  of  R(F)  is  an  interior  point  of  a  bridge  (if  span(x,F)  consists  of  exactly  two 
directions),  an  endpointof  a  bridge  of  positive  length  (if  span(x,F)  consists  of  two 
components  only  one  of  which  has  positive  measure),  or  a  bridge  by  itself  (if 
span  (x,F)  consists  of  two  components  both  of  positive  measure) . 

For  F  labelled  as  in  Lemma  2.3.4  and  H=R(F),  we  define  the  counterclockwise 
traversal  of  dH  to  be  the  circuit  g{f\,h) gifmJi)-  The  counterclockwise  traver- 
sal visits  every  thin  point  twice  and  other  points  once.  Hence  for  x,y  not  thin  points 
we  can  unambiguously  define  dH[x,y]  to  be  the  section  of  the  counterclockwise 
traversal  of  dH  from  x  to  y.  Similarly  if  the  points  in  H'CdH  are  not  thin  points, 
they  have  an  unambiguous  counterclockwise  ordering  on  dH.  If  F  is  extreme,  then 
the  counterclockwise  ordering  covers  every  point  of  F.  Even  if  the  sequence  of 
points  xi,  ...,  xic  on  dH  includes  thin  points,  it  is  still  meaningful  to  say  that  they 
appear  in  counterclockwise  order,  if  they  are  visited  in  that  order  in  a  single  counter- 
clockwise traversal  of  dH.  Of  course,  another  order  differing  in  the  position  of  thin 
points  may  be  consistent  with  a  counterclockwise  traversal  as  well. 

Corollary  2.3.5:  The  center  of  F  is  the  center  of  R  (F) . 

Proof:  It  is  sufficient  to  show  that  for  any  x^U,  rad(x,F)=rad(x,RiF)).   A  point  of 

H=R(F)  maximally  distant  from  x  must  lie  on  dH.    By  the  previous  Lemma  any 

point  on  dH  lies  on  a  geodesic  connecting  two  points  of  F;  hence  by  Lemma  2.2.1, 

the  maximally  distant  point  must  be  an  endpoint  of  the  geodesic,  that  is,  a  point  of  F. 

Thus  a  point  of  H  furthest  from  x  is  necessarily  a  point  of  F  and  rad(x,H)  =  rad(x,F). 

a 

Lemma  2.3.6:  If  /,/'  are  extreme  points  of  F,  then  every  extreme  point  counter- 
clockwise from  /to  /'  is  in  U\f,f'],  and  all  but  /and  /'  lie  in  the  relative  interior  of 
U[f,f']. 

Proof:  The  proof  is  by  induction  on  the  counterclockwise  order  around  F  using  r.  If 
r(f)=^f' ,  then  since  /'  must  lie  in  the  relative  interior  of  U[r(f),f],  r(J)  must  lie  in 
the  relative  interior  of  U[f,f'].  Now  suppose  /"  is  counterclockwise  of /before  /' 
and  r(f")^f' .  By  inductive  assumption,  /"  is  in  the  relative  interior  of  U[f,f']. 
Again  we  have  /'  in  the  relative  interior  of  U[r(f"),f'],  so  r(/")  must  be  in  the 
relative  interior  of  U[f",f'].  Now  U[f",f']  is  contained  in  U\f,f']  except  possibly 
for  a  region  bounded  by  segments  f's'  and  f's"  and  dU{s' ,s"],  where  s'  and  s"  are 
the  closest  shadows  of  g(f,f')  and  g(f" ,f'),  respectively.  But  g(f" ,r(f"))  cannot 
intersect  segment /'j',  since/'  extreme.  Hence  r(f")^U\f,f'].  Since  r(/")  cannot 
lie  on  g(f,f')  as  /,  /',  r(f")  are  extreme  and  distinct,  in  fact  r(/")  is  in  the  relative 
interior  of  f/|/,/'].   □ 

An  immediate  consequence  of  Lemma  2.3.6  is  that  counterclockwise  order  of 
extreme  points  is  an  absolute  order,  not  depending  on  other  extreme  points:  if  a,b,c 
are  extreme  points  of  both  sets  A  and  B,  then  their  order  around  A  is  the  same  as  it  is 
around  B. 

Lemma  2.3.7:  If  points  x,  y,  z,  w  occur  in  that  order  (not  necessarily  all  distinct)  in  a 
counterclockwise  traversal  of  the  boundary  of  a  relatively  convex  polygonal  region  R, 
then  gix,z)  intersects  g{y,w). 
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Proof:  By  creating  dummy  plateaus  if  necessary,  we  can  assume  that  x,  y,  z,  and  w 
each  lie  in  a  plateau.  We  proceed  by  induction  on  the  number  of  plateaus  in  the 
decomposition  of  R.  If  there  is  only  a  single  plateau,  then  /?  is  a  simple  polygon  and 
the  claim  is  immediate.  Otherwise  there  is  some  bridge  g(a,b)  whose  removal  splits 
R  into  relatively  convex  polygonal  regions  Ri  and  Ri-  If  all  of  x,y,2,w  are  contained 
in  one  of  /?i  or  /?2,  the  claim  follows  by  induction.  Otherwise,  since  x,y,z,w  appear 
in  that  order  in  a  counterclockwise  traversal  of  the  boundary  of  R,  there  are  up  to 
symmetry  two  cases:  x,y,a  in  /?i  and  z,w,b  in  R2,  or  x,y,z,a  in  Ri  and  b,w  in  ^2-  In 
the  first  case  gix,z)  and  giy,w)  both  contain  g{a,b).  In  the  second  case  x,y,z,a 
appear  in  that  order  in  a  counterclockwise  traversal  oi  R\,  hence  by  induaion 
hypothesis,  g{x,z)  intersects  giy,a).  Since  giy,a)Cg{y,w),  g(x,z)  must  intersect 
giy.w).   D 

We  say  segment  fxQU  connects  a  relatively  convex  set  i?  to  df/if /is  an  extreme 
point  of  R,  x€dU,  and  fx  intersects  R  only  at  /  (possibly  f=x  or  some  point  of  fx 
besides  x  lies  on  dU). 

Lemma  2.3.8:  (Comiection  Lemma)  Suppose  f\,...,fm  are  extreme  points  of  rela- 
tively convex  set  R  and  segments  fiXi  are  pairwise  disjoint  and  connect  /?  to  dU.  Then 
the  order  of  the  /i's  around  BR  is  the  same  as  the  order  of  x.-'s  around  dU. 
Proof:    Suppose   fi,fj,fk    are    in    counterclockwise    order    around    dR;    we    show 
XjidU[xi,Xfc],  i.e.    Xi,Xj,Xk  are  in  counterclockwise  order  around  dU.    By  Lemma 
2.3.6,  fj  is  in  the  relative  interior  of  C/[/i,/t].   Since  g(Ji,fk)QR  and  fjXj  meets  R  only 
at  fj,  fjXj  cannot  intersect  g(fi,fk)-   Clearly  fj  is  in  the  region  bounded  by  dU[xi,Xk], 
segment  fkXk,  gifkJd .  and  segment  fxi.   Since  Xj^dU and  fjXj  does  not  intersect  fkXk, 
g (fk  Ji)  or  fiXi ,Xj(id U[xi ,Xk]  and  Xj i- x,-  and  Xji'Xk.  a 
Lemma  2.3.9  Mappings  r  and  /  restricted  to  flf/ preserve  order. 
Proof:   Suppose  u\,U2,u-iidU  are  in  that  counterclockwise  order  and  ri  =  r(ui), 
r2=r(M2).  and  r3=r(u3)  are  distinct.  Let  ?,•  be  the  closest  shadow  of  g(M,-,r,). 

We  show  that  Fi,  r^,  Fs  lie  in  that  counterclockwise  order  on  dU\  this  follows 
from  the  claim  that  dU[uj,rj]^dU[Ui,Ti'\  for  all  i^^j.  To  establish  this  claim,  suppose 
for  the  sake  of  contradiction  that  3C/[My,r;]Cat/[u,,r,]  for  some  i#;.  Then 
ry€t/[M,-,r.]  since  u^.r, €£/[«,•, rj,  but  rj(iFCU[ri,Ui]  by  2.3.2(1),  so  rjdjiuiJi).  Jn 
fact  rj^giui.ri)  since  ry  cannot  lie  past  r,-  on  ^(m,-,?,).  If  rj^dU  then  rj=rjidU[ui,ri], 
so  either  ri=rj  or  g(Uj,rj)  can  be  extended  to  giuj.n),  a  contradiction.  If  rjidU, 
then  neither  giui.Ti)  nor  giujjj)  can  bend  at  rj.  Since  Uj,rj€U[ui,ri],  the  two  geo- 
desies must  overlap  in  some  link  containing  r,-.  Furthermore  the  link  must  be 
traversed  in  the  same  direction  in  both  geodesies,  since  Uj,  ui,  r,-,  Tj  appear  in  that 
counterclockwise  order  on  dU.  But  then  g{uj/j)  could  be  extended  to  g(uj,ri),  a  con- 
tradiction. 

We  claim  that  r,?,-  and  rjTj  are  disjoint  whenever  i'^j;  the  Lemma  then  follows 
immediately  by  the  Connection  Lemma.  Suppose  riTi  interserts  tjTj.  The  intersection 
must  be  a  single  point  distinct  from  both  r,-  and  rj.  Hence,  either  rjiU[ri,Ui]  or 
riiU[rj,Uj],  contradicting  2.3.2(1).  o 
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2.4.  Far  Sides 

We  consider  the  notion  of  the  "far  side"  of  a  set  F  from  a  point.  Informally,  an 
extreme  point  f^F  is  on  the  far  side  of  F  from  x  if  the  geodesic  g{x,f)  leaves  R{F) 
after  /.  There  are  two  main  reasons  for  studying  far  sides.  First,  we  use  this  notion 
to  relate  the  order  of  extreme  points  of  F  with  the  order  of  points  on  dU  (Corollary 
2.4.2).  Second,  we  will  be  able  to  compute  relative  convex  hulls  efficiently  using  far 
sides  (Lemma  2.4.3). 

Let  F  be  a  finite  set  of  points  in  U  and  x^dU.  The  far  side  of  F  from  x  is  the  set 
of  all  extreme  points  of  F  counterclockwise  from  r(x)  to  l(x),  inclusive.  Here  r(jc) 
and  l(x)  are  the  clockwise  and  counterclockwise  extreme  points  of  F  from  x,  as 
before. 

Lemma  2.4.1:  Suppose  /is  an  extreme  point  of  F,  x^dU,  and  f^^x.  Then  the  follow- 
ing are  equivalent. 

(1)  /is  on  the  far  side  of  F  from  x. 

(2)  /is  an  extreme  point  of  FU{x}. 

(3)  ff  connects  R(F)  to  dU,  for  any  shadow  /of  g(x,f),  and  there  is  no  geodesic 
extending  g{x,f)  to  an  extreme  point  /'  of  F  distinct  from  /. 

Proof:  If  x^R{F)  then  all  three  conditions  are  satisfied,  so  assume  xiR{F).  Then  x 
is  an  exterior  point  of  F.  Set  r=r{x)  and  l=l{x).  By  Lemma  2.3.2(1) 
f(:U[r,x]f\U[x,l]. 

(1  -  2)  If /is  on  the  far  side  of  F  from  x,  then/€C/[r,/]  by  Lemma  2.3.6.  If /#r  and 
/#/,  then  by  the  Triangle  Lemma,  g(x,f)  intersects  g(l,r)  at  some  point  x' ,  where 
j:'#/ since  /  is  extreme.  Hence  Q(f,x)  =  Q(f,x')^Q(f.gir,l))  C  span(f,F),  and 
spanif  ,FKJ{x})= spanif  ,F) ,  so  /is  extreme  in  FU{ar}.  If /=r,  then  by  the  definition 
of  clockwise  extreme  point,  r{x)  is  unchanged  if  F  is  replaced  with  F\J{x},  so  r  is 
extreme  by  Lemma  2.3.2(3).  The  case/=/  is  similar. 
(2  -  3)  Trivial. 

(3  -  1)  We  argue  the  contrapositive.  If /is  not  on  the  far  side  of  F,  then  by  Lemma 
2.3.6, /€C/[/,r].  Since /cannot  be  on  J(/,r), /€ Arc/.  If  f(:g{x,r)  or /^^(x,/),  then 
g  (.x,f)  can  be  extended  to  an  extreme  point  of  F  distinct  from  /.  Otherwise  by  the 
Triangle  Lemma,  for  some  shadow  /  of  ^(x,/),  gixj)  intersects  g(r,l)  at  a  point  dis- 
tinct from  /,  thus  /f  cannot  connect  R(F)  to  dU.  Q 

Corollary  2.4.2:_Suppose  {/;}  is  a  set  of  distinct  points  on  the  far  side  of  F  from  x 

and,  for  each  i,  /  is  a  shadow  of  g(x,fi).   Then  the  order  of/  on  dU  is  the  same  as 

the  order  of/  on  dR{F). 

Proof:  If  i^j  then  gixji)  and  g(x,fj)  must  diverge  before  reaching  /•  or  fj.    Hence 

fifi  does  not  intersect/}/).  The  result  then  follows  from  the  Connection  Lemma,   a 

Lemma  2.4.3:  Suppose  u,vidU,  u¥=v,  FCU[v,u].  F(Lg{u,v),  and  H=R(F).    Then 

dU[u,v],  g{v,r{v)),  af/[r(v),/(u)],  g(l(u),u)  constitute  a  counterclockwise  traversal 

of  the  boundary  of  R{FUdU[u,v]). 

Proof:  In  this  proof,  we  write  for  example  rn(x)  and  laix)  for  the  clockwise  and 

counterclockwise  extreme  points  of  H  from  x.   This  indicates  the  dependence  upon  H 

explicitly.   Thus  in  the  statement  of  the  Lemma,  r(x)  really  refers  to  rpix).    We  also 
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assume  dU[u,v]^R(F),  hence  some  point  of  dU[u,v]  is  extreme  in  R(FlidU[u,v]). 
If  dU[u,v]CR(F),  a  similar  and  easier  argument  suffices. 

Let  G  be  {m,v}  together  with  the  set  of  convex  comers  of  dU[u,v].  Clearly 
R(F\JdU[u,v])=R{FUG).  We  first  show  that  r/ruG(v)=r/r(v);  we  will  use  Lemma 
2.3.2(5).  Some  point  of  F  is  in  the  relative  interior  of  U[v,u]  by  assumption;  it  is 
easy  to  check  that  in  fact  rfiv)  must  be  in  the  relative  interior  of  U[v,u].  Hence 
FUdU[u,v]CU[rf(v),v].  Any  geodesic  extending  giv.rpiv))  must  stay  in  the  relative 
interior  of  U[v,u],  hence  must  avoid  dU[u,v],  and  must  also  avoid  any  f^F  distinct 
from  rf(v)  since  rf(v)  is  extreme  in  F.  Hence  by  Lemma  2.3.2(5),  r^uG(v)  =  rf(v). 
By  Lemma  2.3.2(3),  rf(v)  is  extreme  in  FUG.  By  a  similar  argument, 
/jruG(")  =  /f(«)  is  extreme  in  FUG. 

We  claim  A=/FuG(''f(v))  is  the  most  clockwise  point  of  G  (ordered  along  dU) 
that  is  extreme  in  FUG:  h  is  either  v,  the  convex  comer  immediately  clockwise  of  v, 
or  possibly  u  if  dU[u,v]  contains  no  convex  comers.  We  have  /i€C/[rf(v),v]  and 
v^U[rF(v),h].  This  is  only  possible  if  v,  h  and  rfiv)  lie  on  a  common  geodesic. 
Since  h  and  rf(v)  are  extreme  in  FUG,  they  must  be  the  endpoints  of  the  geodesic. 
Hence  either  v=/i=/fuG('"F(v)),  or  g(rf{v),v)  extends  to  g(rf(v),h).  Certainly 
A€C/[r/r(v),v].  so  v^dU  and  GCC/[rf(v),v]  force  h^dU[u,vy,  and  g(rF(v),h)  not 
extendible  to  another  point  of  FUG  forces  A  to  be  w  or  the  convex  corner  of  G 
immediately  clockwise  of  v. 

Set/i=rf(v),/2=rf(/i),  ....  fh=lF(u)=rF(fh-i).    Let  G'  =  {gi gi}  be  the 

subset  of  G—F  extreme  in  FUG,  where  the  index  is  given  by  the  ordering  along 
dU[u,v].  We  want  to  show  that  r/ryo  maps  /i  -/i  -  •  •  •  -/a  -  ^i  -»  •  •  •  ■*  gi -"fi- 
We  have  already  shown  that  lf\jG  maps  f\  to  gi,  hence  rF\jG(gi)=fi.  By  a  similar 
argument  rfuG(//.)=5i-  It  is  easy  to  see  that  rFuaigd  can  only  be  ^,+i,  for  l<j</. 

It  remains  to  show  that  rF[jG(fi)=fi+i,  for  l<i<h.  This  follows  if  we  establish 
span(fi,FUG)=span(fi,F)  (which  also  establishes  fi  extreme  in  FUG)).  Since 
fh=lF(u)  is  extreme  in  FUG,  /^  is  extreme  in  FU{v},  hence  fn  is  on  the  far  side  of  F 
from  V,  and  hence  f\,....,fh  are  on  the  far  side  of  F  from  v.  Similarly  f\,...,fh  are  on 
the  far  side  of  F  from  u.  Since  r  and  /  preserve  order  on  dU,  the  far  side  of  F  from 
any  w^dU[u,v]  includes  f\,...,fh-  In  particular,  for  any  w^dU[u,v]  and  any  i, 
\<i<h,  we  have  /i«g(/(w),r(w))  (since  /)•  extreme  in  F  and  fii'l{w),r{w)),  and  (as 
in  the  proof  of  Lemma  2.4.1),  8(y;-.w)€e(/i.5(/(w),r(w)))  C  Q(fi.R(F))  = 
span(fi,F).   Hence  span(fi,F\JG)  =  span(Ji,F). 

Now  if  gi  =  u  and  gi=v,  we  are  done.  If,  say,  gi^u,  then  u  must  appear  on 
g(fk,gi)-  To  obtain  the  Lemma  split  g(fk,gi)=gilF(u),gi)  into  gihM.u)  and 
^(m.^i),  then  merge  ^(u.^i)  with  dU[gi,gk].   A  similar  split  applies  if  g/#v.   □ 

2.5.  The  General  Position  Assumption 

For  the  remainder  of  the  paper  we  make  the  following  general  position  assump- 
tion: no  reflex  corner  of  dU  is  equidistant  from  two  sites.  This  condition  can  always 
be  satisfied  by  applying  a  slight  perturbation  to  the  positions  of  the  sites  or  corners. 
If  this  assumption  is  not  made,  the  set  of  points  equidistant  from  two  sites  may 
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include  a  two-dimensional  region  [A87].  This  seaion  contains  some  consequences  of 
the  general  position  assumption  that  are  critical  to  the  rest  of  the  paper.  We 
emphasize  that  none  of  Lemmas  2.5.1,  2.5.2,  and  2.5.3  hold  if  the  assumption  is 
removed. 

Lemma  2.5.1:  If  s  and  t  are  distinct  sites,  x^U,  and  d,(x)=d,(x),  then  Q(x,s)=^B(x,t). 
Proof:  Suppose  Q(x,s)  =  Bix,t).  Then  the  geodesies  g(s,x)  and  g(t,x)  share  their  final 
link  yx.   Point  y  must  be  the  anchor  of  x  with  respect  to  both  of  s  and  t.   Now  y^s,t 
since  otherwise  d(s,y)=d(t,y)  would  imply  s=t.   Hence  y  must  be  a  reflex  comer  of 
at/ equidistant  from  s  and  t,  contradicting  the  general  position  assumption,  n 

Lemma  2.5.2:  Suppose  u,v€U,  u^v,  and  each  of  u  and  v  is  equidistant  from  distinct 
sites  s  and  t.  Then  g(s,u)  does  not  intersect  g(t,v). 

Proof:  Suppose  x  lies  on  both  gis,u)  and  ^(r,v).  Without  loss  of  generality,  assume 
d(x,s)^d(x,t).  Observe  that 

dy(t)=dy(s)^dy(x)+d(x,s)^dy(x)+d(x,t)=dyit). 

Hence  d{x,s)=d{x,t),  path  g{s,x)KJg{x,v)  is  in  fact  geodesic  g{s,v),  and  similarly 
g(,t,x)\Jg{x,u)  is  g(t,u).  Since  x  lies  on  both  geodesies  g{s,u)  and  g(s,v),  x  must  be 
a  reflex  comer  of  3  C/  or  an  interior  point  of  a  common  link  of  the  two  geodesies  con- 
necting two  reflex  comers  of  dU.  In  the  second  case  the  link  must  be  common  to 
g(t,u)  and  g(t,v)  as  well,  hence  in  either  case  we  can  find  a  reflex  corner  of  dU 
equidistant  from  s  and  t,  violating  the  general  position  assumption.  □ 

Lemma  2.5.3:  There  is  at  most  one  point  equidistant  from  three  distinct  sites. 
Proof:  Suppose  to  the  contrary  that  points  u  and  v  are  both  equidistant  from  sites 
r,s,t.   First  note  that  r,s,t  cannot  lie  on  a  common  geodesic,  for  if  say  r€g(s,t)  then 
by  Lemma  2.2.1,  du{r)<max(du(s) ,du{t)) .   Hence  r,s,t  are  extreme  points  of  {r,s,t}. 

We  claim  that  r,s,t  are  extreme  points  of  {r,s,t,u}  (and  also  of  {r,s,t,v}).  In 
order  to  demonstrate  this,  we  show  riR({s,t,u}).  Now  r  does  not  lie  on  g(u,s)  or 
g{u,t),  else  5  or  r  would  be  further  from  u  than  r.  As  argued  before,  r  does  not  lie 
on  g(.s,t).  If  r  is  in  the  interior  of  R({s,t,u})  then  by  the  Triangle  Lemma,  g{u,7) 
intersects  g(s,t),  where  ?  is  a  shadow  of  g{u,r).  But  then  using  Lemma  2.2.1  again, 
we  would  have  rf„(r)<max(<i„(j),J„(r)). 

Now  suppose  one  of  u  and  v,  say  u,  is  extreme  in  R{{r,s,t,u,v}).  Hence  r,s,t,u 
are  extreme  in  R{{r,s,t,u))\  assume  that  they  appear  in  that  counterclockwise  order. 
By  Lemma  2.3.6,  r^U[u,s\,  r€£/[j,u],  and  r,tigi,u,s).  It  must  be  that 
giu,s)  =  giu* ,7)  intersects  either  giv,r)  or  g(v,t),  assume  it  is  g(v,t).  To  obtain  a 
contradiction  of  Lemma  2.5.2,  we  show  in  fact  g{u,s)  intersects  ^(v,:).  Now  uu* 
intersects  R({r,s,t,u,v})  only  at  u  since  u  is  extreme  in  R({r,s,t,u,v}).  Hence  uu*  is 
disjoint  from  /?({r,J,^v}).  Also  jJ intersects  R({r,s,t,v})  only  at  s,  since  s  is  extreme 
in  Ri{r,s,t,v})  and  some  portion  of  g{u,s)  must  lie  in  R{{r,s,t}),  hence  in 
R({r,s,t,v})-   Since  g(v,t)CR({r,s,t,v}),  giv,t)  must  intersects  g(u,s). 

If  neither  u  nor  v  is  extremal  in  R({r,s,t,u,v})  then  both  u,viR({r,s,t})  and  the 
proof  is  similar.  Geodesies  g(u,r),  g(u,s),  and  g(u,t)  split  Ri{r,s,t})  into  three  geo- 
desic  triangles.     Hence    v   lies    in    one    of   the    triangles    which    again    implies    a 
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contradiction  of  Lemma  2.5.2.   □ 

2.6.  Voronoi  CelTs 

The  bisector  b{s,t)  of  distinct  sites  s  and  /  is  {xiU:d,(x)=d,(x)}  and  the  half 
space  closer  to  s,  H(s,t),  is  {x€U:ds(x)<d,ix)}.  Clearly  H(s,t),  b(s,t),  and  H(t,s) 
form  a  partition  of  U.  A  breakpoint  of  bis,t)  is  the  intersection  of  b(s,t)  with  a  shor- 
test path  partition  edge  from  s  or  t.  Figure  2.6.1  indicates  the  bisectors  of  three 
points,  with  breakpoints  marked. 


Figure  2.6.1.  Bisectors;  dots  indicate  breakpoints 

Lemma  2.6.1:  [A87,  3.22]:  Bisector  b(s,t)  is  a  smooth  curve  connecting  two  points 
on  dU  and  having  no  other  points  in  common  with  dU.  It  is  the  concatenation  of 
0(n)  straight  and  hyperbolic  arcs;  the  points  where  such  arcs  meet  are  precisely  the 
breakpoints  of  bis,t).  The  tangent  to  b{s,t)  at  x  bisects  the  angle  between  B(x,s)  and 
Bix,t). 

In  particular,  together  with  Lemma  2.5.1  this  implies  that,  given  x^b(s,t), 
Qix,s)  and  -9(;c,r)  enter  H(s,t)  at  x,  while  -9(.x,s)  and  Q(x,t)  enter  H(t,s)  (if  they 
stay  within  U,  that  is). 

Corollary  2.6.2:  The  relative  boundary  of  H(s,t)  and  H{t,s)  is  bis,t). 

Lemma  2.6.3:  [A87,  3.17]  H{s,t)  is  path-connected. 

Recall  that  S  is  the  set  of  sites.  The  (geodesic  furthest-site)  Voronoi  cell  of  site  s 
is  Vis)=  nHit,s).  The  (geodesic  furthest-site)  Voronoi  diagram  V  is  {x€b(s,t):s,t^S 

and  dj(x)=max  dr(x)}^   Figure  2.6.2  indicates  the  Voronoi  diagram  and  Voronoi  cells 

of  the  three  points  depicted  in  Figure  2.6.1.  A  Voronoi  edge  e{s,t)  is  Vr\b(s,t),  if 
the  intersertion  consists  of  more  than  one  point  (else  we  will  say  that  e(,s,t)  does  not 
exist) .  A  (Voronoi)  vertex  is  a  point  x^V  which  has  three  or  more  sites  furthest  from 
it.  By  Lemma  2.5.3  above,  there  is  at  most  one  such  point  x  for  each  triple  of  sites. 
A  hitpoint  is  the  intersection  of  a  Voronoi  edge  with  dU.  Intuitively,  a  hitpoint 
correspond  to  the  "point  at  infinity"  of  an  infinite  Voronoi  edge  in  a  Euclidean 
furthest-site  Voronoi  diagram. 

Lemma  2.6.4:  Voronoi  edge  e(^s,t)  is  connerted  and  has  vertices  or  hitpoints  as  end- 
points. 
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Figure  2.6.2.  Voronoi  cells 
Dashed  lines  appear  in  refined  Voronoi  diagram 

Proof:  Suppose  r,s,t  are  distinct  sites.  Since  dr  is  continuous,  a  connected  component 
of  b(s,t)r\(b(r,s)UH(r,s))  must  have  for  each  of  its  endpoints  either  an  endpoint  of 
b(s,t),  i.e.,  a  hitpoint,  or  a  point  equidistant  from  s,t,r.  But  by  Lemma  2.5.3,  there 
is  at  most  one  point  equidistant  from  s,t,r,  so  b(s,t)n(b(r,s)UH(r,s))  consists  of  a 
single  connected  component.   Hence  also  eis,t)=b{s,t)n   H  (bir,s)UH(r,s))  is  con- 

r*s,t 

nected  and  has  hitpoints  or  vertices  for  endpoints.   o 

Lemma  2.6.5:  Suppose  si,...,sic,sic+i=si  are  the  sites  furthest  (and  thus  equidistant) 
from  vertex  v,  and  directions  9(v,ji),...,9(v,5t)  are  in  counterclockwise  order.  Then, 
for  each  i,  edge  e(Si,Si+i)  is  incident  to  v  and  extends  away  from  v  in  direction 
bisecting  Z.J,+iV5,-,  as  long  as  that  direction  (locally)  stays  inside  U. 
Proof:  Elementary  analysis,  using  V(i,.(v)=  —  9(v,j,).  a 

If  vertex  v  appears  on  dU,  then  there  is  only  one  edge  of  V  incident  to  v:  as  v 
cannot  be  a  comer  of  dUhy  the  general  position  assumption,  it  must  be  an  interior 
point  of  a  wall.  Hence  only  the  edge  bisecting  Lsvt  remains  within  U,  where  direc- 
tions 9(v,j)  and  9(v,r)  are  the  most  clockwise  and  most  counterclockwise  directions 
towards  sites  furthest  from  v,  respectively. 

Lemma  2.6.6  (Extension  Lemma):  If  x  lies  on  gis,y)  and  x^V(s)  or  x^eisj),  for 
some  other  site  t,  then  all  of  g(s,y)  past  x  lies  in  V{s). 

Proof:  Suppose  x^V(s).  For  any  site  r^^j  and  for  any  z^g(x,y),  d^iz)  = 
d,(x)+d(x,z)  >  dr(x)+d(x,z)  2:  driz),  so  z^V{s).  Suppose  x^e(s,t)  and  y=^x.  By 
the  general  position  assumption  x  is  not  a  comer  so  Qix,y)=  —Q(x,s).  Since  y^'x, 
direction  Qix,y)  stays  in  t/  at  x.  We  show  d{x,y)  enters  V(s)=  DHir.s)  at  x;  the 

result  follows  as  before.  Since  x^e(s,t),  for  any  r#j,  dr(x)^dj(x).  If  dr(x)=ds(x), 
-Q{x,s)  =  Qix,y)  enters  Hir.s);  if  drix)< d six),  any  direction  (locally)  stays  in  H{r,s). 
a 

An  immediate  consequence  of  this  Lemma  is  that  every  point  in  a  Voronoi  cell 
is  connected  to  dU:  if  x^V(s),  then  segment  xxCV(s),  where  x  is  a  shadow  of  g(s,x). 

Lemma  2.6.7:  Both  V(s)r\dU  And  V(s)  are  path-connected. 

Proof:  Since  every  point  of  V(s)  is  connected  to  a  point  of  V(s)ndU,  it  suffices  to 
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show  V(s)r\dU  is  connected.  To  prove  it,  we  will  show  that  H{r,s)r\H(t,s)ndU  is 
connected  for  every  r,t¥'s.  Label  the  hitpoints  of  b(r,s)  and  b(t,s)  as  Xr,yr  and  x,,y,, 
respectively,  so  that  (dUnH(r,s))\J{xr,yr}  =  dU[xr,yr]  and  (dUnH(t,s))U{x„y,}  = 
dU[x,,y,].  The  only  counterclockwise  ordering  of  these  points  that  disconnects 
dUnHir,s)nH(t,s)  is  Xr,  yt,  x,,  yr-  In  particular,  this  implies  that  H{r,s)\JH{t,s) 
cover  dU.  We  show  this  is  impossible;  suppose  this  were  the  ordering.  Since 
siH(,r,s)UH(t,s),  Hir,s)UH(t,s)=^U,  so  bir,s)  must  intersect  b{t,s)  in  at  least  two 
distinct  points.  But  this  contradicts  Lemma  2.5.3,  since  each  of  the  two  points  would 
be  equidistant  from  sites  r,  s,  and  t.   d 

2.7.  The  Ordering  Lemma 

Let  C  be  the  relative  convex  hull  of  S,  the  set  of  sites.  The  first  Lemma  of  this 
section  shows  that  if  V{s)  is  not  empty,  then  s  must  be  an  extreme  point  of  S  (and  in 
fact  on  the  far  side  of  S  from  any  point  in  the  closure  of  V(s)ndU).  Hence  we  can 
assume  that  all  sites  are  extreme  and  that  they  are  ordered  by  the  counterclockwise 
traversal  of  dC.  The  main  result  of  this  section  is  the  Ordering  Lemma:  the  order  of 
Voronoi  cells  around  dU  is  exactly  the  order  of  sites  around  dC.  In  addition  we  show 
that  there  is  a  collection  of  at  most  three  geodesies  separating  every  point  of  dU  from 
its  furthest  site  (or  sites).  This  is  used  to  prove  the  rather  remarkable  fact  that  there 
are  at  most  0(n+k)  distinct  links  in  all  geodesies  that  connect  corners  of  dU  to  their 
respective  furthest  sites.  We  remark  that  Lemmas  2.7.2,  2.7.5,  and  2.7.6  are  based 
on  very  similar  Lemmas  proved  by  Suri  [S87]. 

Lemma  2.7.1:  If  V{s)  is  not  empty,  then  j  is  an  extreme  point  of  S  on  the  far  side  of 
C  from  any  x  in  the  closure  of  V(s)ndU. 

Proof:  Suppose  5  is  a  site  furthest  from  x.  If  s€gir,t),  for  r,t  sites  distinct  from  s, 
then  by  Lemma  2.2.1,  d,{x)<m.aj(.(dt(x) ,dr(x)) ,  a  contradiction.  Suppose  the  last  seg- 
ment of  g{x,s)  can  be  extended  beyond  s  staying  in  C;  let  s'  be  a  point  lying  on  such 
an  extension  and  on  the  boundary  of  C.  Then  d(s' ,x)>dj(x),  s'  must  lie  on  g(t,u) 
for  some  sites  t,u,  and  d(s' ,x)^rmixid,(x),duix)),  contradicting  the  choice  of  s. 
Hence  s  is  not  an  interior  point  of  C,  so  s  must  be  an  extreme  point.  By  Lemma 
2.4.1,  s  is  on  the  far  side  of  C  from  x.   n 

Lemma  2.7.2:  Suppose  s  and  t  are  furthest  sites  from  u,v^U,  respectively,  with  m^v, 
s^t.  Then  g{u,t)  does  not  meet  g(v,s). 

Proof:  Suppose  to  the  contrary  that  x^g(u,t)ng(v,s).  By  the  triangle  inequality  we 
have  d(s,u)^d(s,x)+d(x,u)  and  d(t,v)^d(t,x)+d(x,v).  Adding  we  obtain 
d(s,u)+d(t,v)<d(t,u)+d(s,v).  However,  since  s  and  t  are  furthest  sites  from  u  and 
V  respectively,  d(s,u)^d(t,u)  and  dit,v)^d(s,v).  Thus  we  must  have  d(t,v)  =  d(s,v) 
and  d(s,u)=d(t,u).  However,  in  this  situation  g(u,t)  cannot  intersect  ^(v,^)  by 
Lemma  2.5.2,  a  contradiction,  o 

Lemma  2.7.3:  Suppose  V(r) n 3 C/ immediately  follows  V(s)ndU  in  counterclockwise 
order  around  dU.   If  ui^sj  is  another  site  on  dC[s,t],  then  V(«)  is  empty. 
Proof:  Let  x  be  the  hitpoint  of  e{s,t)  so  that  (near  jc)  V{s),  x,  V{t)  are  in  counter- 
clockwise order  along  dU.    Let  7  and  t  be  closest  shadows  of  g{x,s)  and  g{x,t). 
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respectively.  We  show  that  x,  J,  t  are  distinct  and  appear  in  that  order  ^counterclock- 
wise around  dU.  Note  that  by  Lemma  2.5.1,  g{x,'s)=g{x,s)  and  ^(j:,r)=J(x,r)  are 
geodesies  emanating  from  x  with  distinct  initial  direaions;  hence  7#r.  Near  x, 
Vit)ndU=H(s,t)ndU  and  V(s)naU=Hit,s)ndU;  also  Q(x,s)  enters  H{s,t)  and 
Qix,t)  enters  H(t,s).  Since  g(x,s)  does  not  intersect  gix,t)  again,  the  ordering  of  x, 
J,  t  must  be  counterclockwise  around  dU. 

We  claim  7,  u,  and  r  are  in  that  counterclockwise  order  on  dU,  where  u  is  the 
closest  shadow  of  g{x,u).  Let  r  and  /  be  the  clockwise  and  counterclockwise  extreme 
points  of  5  from  x,  with  r  and  7 closest  shadows  of  g{x,r)  and  ^(x,/),  respectively. 
By  Lemma  2.3.2(2),  aC/[r,7]  does  not  contain  x.  By  Lemma  2.7.1,  s,t  are  on  the  far 
side  of  S  from  x.  Hence  by  Corollary  2.4.2,  7,r€aC/l7,7],  and  in  fact  the  counter- 
clockwise order  must  be  r  ,  7,  r,  7,  since  aC7[7,r]  does  not  contain  x.  By  Corollary 
2.4.2,  r,  s,  t,  I  appear  in  that  counterclockwise  order  on  C;  since  u  appears_between  s 
and  r,  u  is  on  the  far  side  of  5  from  x.  Again  by  Corollary  2.4.2,  7,  u,  t  appear  in 
that  counterclockwise  order  on  dU. 

Now  suppose  contrary  to  the  Lemma  V(m)  is  not  empty;  we  will  obtain  a  con- 
tradiction. Since  Voronoi  cells  are  connected  to  dU,  there  is  y€V(u)n3C/.  Now  y#x 
since  Voronoi  cells  are  relatively  open,  and  x  lies  on  the  boundary  of  V(,s)  (and 
V(r)). 

As  s,  u,  and  t  are  distinct  extreme  elements  of  S,  the  set  {s,u,t}  is  extreme.  As 
x^dU,  it  is  impossible  for  x  to  lie  in  the  interior  of  ^sut.  Moreover,  s,  u,  and  r  are 
on  the  far  side  of  {s,u,t}  from  x  (in  that  counterclockwise  order).  Thus  either 
xig{s,t)  or  {s,u,t,x}  is  extreme  (with  s,  u,  t,  x,  occurring  in  that  counterclockwise 
order).  In  either  case  g{x,u)  intersects  g{s,t)  and  enters  t/[j,r]  at  some  point  u'  (as 
M  is  an  extreme  point  of  C  lying  counterclockwise  between  s  and  r).  Assume 
y(iU[x,u],  the  case  yiU[u,x]  is  similar.  Lemma  2.7.2  implies  that  g{y,t)  does  not 
intersect  g{x,u).  As  t  lies  in  the  relative  interior  of  U{u,x]  and  y  lies  in  C/[x,m], 
g{y,t)  must  intersect  mu  at  some  point  u"^u.  This  implies  that  the  portion  of  g{x,u) 
from  u'  to  «"  inclusive  is  contained  in  the  triangle  ^yst;  in  particular  u^Ayst.  By 
Lemma  2.2.1,  this  contradicts  the  choice  of  m#j,  r  as  a  site  furthest  from  y.   n 

Corollary  2.7.4  (Ordering  Lemma):  The  ordering  of  sites  with  nonempty  Voronoi 
cells  around  dC  is  the  same  as  the  ordering  of  Voronoi  cells  around  dU. 

Lemma  2.7.5:  Suppose  {ui,...,u„}  and  {vi,...,v„}  are  separated  by  boundary  geo- 
desic gix,y).  Then  there  is  a  total  of  at  most  0(m+n)  distinct  links  in  paths  ^(u,-,v,), 
i  =  l,...,m. 

Proof:  By  a  proof  essentially  identical  to  Lemma  4  of  Suri  [S87],  for  each  i  there  are 
at  most  three  links  of  g(«,-,v,)  not  in  r(x)UT(y).  (Three  links  are  needed  rather  than 
one  because  u,v  need  not  be  comers  of  dU.)  n 

Lemma  2.7.6:  There  are  three  boundary  geodesies  gi,g2fg3  so  that  for  any  point 
a^dU  and  any  site  b  furthest  from  a  one  of  the  geodesies  separates  a  from  b. 
Proof:  Pick  x€5  arbitrarily.   Let  y  be  a  site  furthest  from  x  and  z  a  site  furthest  from 
y.   We  argue  the  case  that  x,y,z  are  distinct  and  appear  in  that  counterclockwise  order 
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around  dC;  the  case  that  x  =z  is  easier  and  the  case  that  the  order  is  x,z,y  is  similar. 

We  claim  that  we  can  choose  a  site  w  furthest  from  z  in  dC[x,y].  We  first  show 
that  w  can  be  chosen  to  lie  in  dC[x,z]:  if  we  can  only  choose  a  w  furthest  from  z  in 
dC[z,x]-{x},  then  in  particular  diz,w)>d(z,y).  We  also  have  d(x,y)^d(x,z)  and 
d(y,z)^diy,w),  so  adding  all  three  we  get  d(x,y)+d(z,w)>d(x,z)+d(y,w).  This 
contradias  the  triangle  inequality  as  can  be  seen  by  considering  a  point  in  the  inter- 
section of  giy.w)  and  g{x,z).  Now  w  cannot  be  in  dC[y,z]-{y},  else  g(x,w)  inter- 
sects g(y,z),  a  contradiction  of  Lemma  2.7.2,  taking  u=x,  t=w,  v=z,  and  s=y. 
Hence  we  can  choose  w  in  dC[x,y]. 

Let  X* ,y* ,z*  and  y,J,w  be  the  closest  foreshadows  and  closest  shadows  of 
gix,y),  giy.z),  and  giz,w),  respectively.  We  claim  y* idU[x* ,y]:  this  is  immediate  if 
y^dU  or  gix,y)  and  g(z,y)  share  a  common  final  link;  otherwise  it  follows  from 
Lemma  2.2.4(1),  since  y*  is  the  closest  shadow  of  g(z,y).  Similarly  z*^dU{y* ,1]  and 
x'^dU[z\w]. 

Now  by  the  Extension  Lemma,  jc*,  y* ,  z*  have  y,z,w  as  their  respective  furthest 
sites.  By  the  Ordering  Lemma,  all  points  in  dU[x* ,y*]  have  furthest  sites  in  dC\y,z]; 
clearly  g{x,y)  separates  dU[x* ,y*]  from  dC\y,z].  Similarly  the  sites  furthest  from 
dU\y' ,z*'\  lie  in  dC[z,w\  and  these  two  sets  are  separated  by  g(,y,z).  Finally,  the 
sites  furthest  from  dU[z* ,x']  lie  in  dC[w,y]  and  these  two  sets  are  separated  by 
J(z,x).  Note  that  in  fact  x*  can  have  furthest  sites  both  in  dC\y,z]  and  BC[w,y],  and 
similarly  for  y*  and  z',  while  each  of  the  remaining  points  oi  dU  has  its  furthest  sites 
in  exactly  one  of  dC\y,z],  dC{z,w],  dC[yv,y].  □ 

2.8.   The  Refined  Voronoi  Diagram 

The  Voronoi  diagram  V  clearly  has  0{k)  edges,  since  together  with  dU  li  forms 
a  planar  graph  with  at  most  k  bounded  regions  and  all  Voronoi  vertices  have  degree 
three  or  more.  However,  this  is  not  an  accurate  description  of  the  size  complexity  of 
the  Voronoi  diagram,  since  each  Voronoi  edge  may  consist  of  sections  of  several  dif- 
ferent hyperbolic  arcs.  This  section  discusses  a  refinement  of  the  Voronoi  partition, 
obtained  by  further  subdividing  each  Voronoi  cell  V{s)  by  the  shortest  path  partition 
from  s.  Each  bounding  edge  of  a  refined  Voronoi  cell  is  a  line  segment  or  a  section 
of  a  single  hyperbolic  arc.  The  main  theorem  is  a  linear  bound  on  the  size  complex- 
ity of  the  refined  Voronoi  diagram.  This  implies  an  0{n+k)  bound  on  the  size  com- 
plexity of  the  Voronoi  diagram  itself. 

The  refined  Voronoi  cell  of  site  s  with  anchor  a,  V'a(s),  is  V(s)r\Pais).  The 
refined  bisector  edge  e'ab(s,t)  is  eis,t)r\Pais)r]Pbit).  The  refined  partition  edge  (fi-om 
s  with  anchor  a),  plis),  is  V(s)npa(s).  A  refined  Voronoi  edge  is  a  refined  bisector 
edge  or  a  refined  partition  edge.  Observe  that  distinct  refined  Voronoi  edges  are  dis- 
joint (except  possibly  at  their  endpoints).  Empty  refined  Voronoi  cells  and  edges  and 
refined  Voronoi  edges  consisting  of  a  single  point  are  disregarded. 

Suppose  e*ai,(s>t)  is  not  empty.  It  is  easy  to  see  that  each  endpoint  of  eabi^'O  is 
either  a  vertex  of  V,  a  hitpoint,  or  a  breakpoint.   Moreover,  e*ab(s,t)  does  not  contain 
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breakpoints  (except  possibly  as  endpoints).   Consequently  g^(5,r)  is  a  hyperbolic  arc 

or  a  line  segment. 

Lemma  2.8.1:  Either  p*a(s)  is  empty,  or  it  is  all  of  pais),  or  it  has  an  open  endpoint 

at  a  breakpoint  of  e(s,t)  for  some  site  t  and  closed  endpoint  on  dU.   (The  latter  two 

cases  are  illustrated  in  Figure  2.6.2.) 

Proof:  Suppose  plis)  is  not  empty  and  is  not  all  of  pais)= ay  (where  y^dU).    Then 

Pais)  must  intersect  some  edge  eis,t)  before  first  entering  Vis).    By  the  Extension 

Lemma,  the  intersection  is  a  single  point  x,  and  p*ais)=xy-{x}  is  contained  in  Vis). 

a 

Suppose  p'ais)  is  a  refined  partition  edge  of  Vis).  Let  tiplis))  be  the  first  link 
of  the  geodesic  gia,s),  directed  towards  s.  Then  by  definition  of  the  refined  parti- 
tion edge  p'ais).  Pais)  and  tiplis))  are  coUinear  and  meet  at  a.  We  claim  that  if 
p'ais)  and  p'a-is')  are  distinct  refined  partition  edges,  then  tiplis))  and  tip'a'is'))  are 
distinct,  at  least  as  directed  links.  This  claim  certainly  holds  if  a#a'.  If  a=a' ,  then 
we  must  have  s=^s'  by  the  previous  Lemma.  Since  p'ais)CVis)  and  Pa'is')CVis'), 
we  must  have  Qia,s)¥'Qia,s').  Since  tiplis))  and  tiplis'))  have  these  two  direc- 
tions, they  must  be  distinct. 

Lemma  2.8.2:  Link  tiplis))  is  a  link  of  giv,s).  where  v  can  be  chosen  to  be  a  corner 
lying  in  dUOVis)  or  a  hitpoint  of  Vis)  (i.e.  an  endpoint  of  dUDVis)). 
Proof:  Let  pais)=ay,  y^dUDVis).  Segment  ay  partitions  U  into  two  polygonal 
regions  Ui  and  U2  so  that  Ui  contains  s  and  every  geodesic  to  s  from  a  point  in  U2 
contains  tiplis)).  If  we  traverse  dUCMJi  starting  at  y,  we  must  encounter  either  the 
endpoint  of  dUC\Vis)  or  a  comer  of  dU  lying  in  Vis).  Q 

The  refined  Voronoi  diagram,  V' ,  is  the  union  of  all  refined  Voronoi  edges.  A 
vertex  of  V*  is  a  vertex  of  V  or  a  breakpoint.  A  hitpoint  of  V'  is  a  hitpoint  of  V  or 
the  point  of  intersection  of  a  refined  partition  edge  with  3C/  (in  the  case  when  the 
refined  partition  edge  coincides  with  a  shortest  path  partition  edge,  only  its  non- 
anchor  endpoint  is  considered  a  hitpoint). 

Lemma  2.8.3:  There  are  at  most  0(n  +Jt)  refined  Voronoi  edges  and  vertices. 
Proof:  Clearly  V  has  k  cells,  hence  k  hitpoints  and  Oik)  edges.    We  show  that  there 
are  only  Oin-¥k)  refined  partition  edges.    Since  each  refined  partition  edge  contri- 
butes a  single  breakpoint,  it  follows  that  there  are  only  Oin+k)  refined  bisector 
edges.  By  planarity,  there  are  only  0(n  +Jt)  vertices  in  V'  as  well. 

By  Lemma  2.7.6,  there  are  three  boundary  geodesies  ^i,  gi,  g-i  and  a  partition 
of  dU  into  three  fragments  U\,  Uj,  U-i  so  that  for  1  =  1,2,3,  gi  separates  every  point 
in  Ui  from  its  furthest  sites.  (This  is  not  strictly  true  for  the  endpoints  of  Ui,  but  the 
argument  is  similar.)  Let  W,-  be  the  union  of  the  sets  of  links  of  geodesies  g(v,w), 
where  v  is  a  corner  of  dU  lying  in  Ui  and  w  is  the  unique  site  furthest  from  v  or 
where  v  is  a  hitpoint  of  V  lying  in  f/,-  and  w  is  one  of  the  two  sites  whose  Voronoi  cell 
boundary  contains  v.  Since  there  are  k  hitpoints  and  n  comers  in  all  of  dU,  there  are 
certainly  at  most  as  many  in  each  C/,-.  By  Lemma  2.7.5,  there  are  Oin+k)  links  in 
Wf.  Now  modify  W,-  so  that  it  contains  for  each  link  two  oppositely  directed  links; 
this  doubles  its  size.    By  Lemma  2.8.2,  if  Pa  is  a  refined  partition  edge,  then  t(pa)  is 
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in  Wi,  for  some  j  =  1,2,3.  Since  each  t(pa)  corresponds  to  a  unique  refined  partition 
edge,  there  are  at  most  0(n  +k)  such  edges.  □ 

2.9.  Directing  Edges  of  tlie  Refined  Voronoi  Diagram 

We  let  c  be  the  center  of  C.  Necessarily  c  lies  on  some  Voronoi  edge,  since 
there  must  be  (at  least)  two  sites  attaining  the  maximum  distance  from  c. 

Lemma  2.9.1:  Suppose  5  is  a  site  furthest  from  x^U.  Then  the  angle  between  Q(x,s) 
and  9(x,c)  is  less  than  it/2. 

Proof:  Suppose  the  angle  between  B(x,s)  and  9(x,c)  is  at  least  11/2.  By  Lemma 
2.2.2,  d(c,s)>dix,s).  But  since  5  is  a  furthest  site  from  x,  rad(x,C)=d{x,s).  This 
contradicts  the  choice  of  c  as  the  center  of  C.  n 

We  use  this  Lemma  to  direct  each  refined  Voronoi  edge  towards  c  (if  there  is  a 
Voronoi  edge  containing  c  in  its  interior,  we  split  the  edge  at  c).  We  show  how  to 
direct  an  edge  e(s,t)  of  the  (unrefined)  Voronoi  diagram;  this  direction  extends  to 
each  refined  bisector  edge.  A  similar  argument  directs  each  refined  partition  edge. 
Suppose  xieis,t),  x#c.  Since  e(s,t)  bisects  Lsxt  and  both  the  angle  between  6(x,c) 
and  B{x,s)  and  the  angle  between  6(x,c)  and  B(x,t)  are  less  than  ir/2,  it  is  not  possi- 
ble that  6(x,c)  is  perpendicular  to  e(s,t)  at  x.  Hence  we  can  direct  eis,t)  towards  c 
locally  at  x.  If  cie{s,t),  this  direction  extends  globally  to  e(s,t)  (since  direction 
towards  c  is  a  continuous  function  away  from  the  comers  and  can  never  be  perpen- 
dicular to  e{s,t)).  If  c^e(s,t),  then  eis,t)  is  split  at  c  into  two  pieces,  each  of  which 
is  consistently  directed  towards  c. 

Lemma  2.9.2:  Suppose  vidU  is  a  vertex  of  V* .  Then  there  are  at  least  two  edges  of 
V*  entering  v  and  exactly  one  edge  of  V*  leaving  v. 

Proof:  We  consider  the  case  that  v  is  in  fact  a  vertex  of  V;  the  case  that  v  is  a  break- 
point is  similar.  For  simplicity  assume  there  are  exactly  three  sites  equidistant  from 
V.  Label  the  sites  r,s,t  so  that  9(v,r),  B(v,j),  and  9(v,r)  are  in  counterclockwise 
order  leaving  v  so  that  mLrvt<'u.  (This  is  possible  by  the  previous  Lemma  as  each  of 
the  angles  formed  between  9(v,c)  and  9(v,r),  9(v,j)  or  9(v,r)  has  measure  less  than 
Kll.)  By  Lemma  2.6.5,  Voronoi  edges  e(r,s),  e{s,t),  and  e{t,r)  are  incident  to  v 
and  extend  in  directions  bisecting  angles  Lsvr,  Ltvs,  and  Lrvt,  respectively.  Hence 
edges  e{r,s)  and  e(s,t)  enter  v  and  e(t,r)  leaves  v.  The  general  case  of  a  vertex  of 
arbitrary  high  degree  is  handled  analogously.  Vertices  vidU  of  degree  two  or  less 
do  not  occur  by  definition  of  V*.  □ 

Suppose  Voronoi  edge  e(s,t)  intersects  dU  at  hitpoint  x;  recall  that  by  the  gen- 
eral position  assumption  x  is  not  a  comer  of  dU.  Then  e{s,t)  is  directed  into  the  inte- 
rior of  U,  because  9(x,c)  makes  an  angle  strictly  less  than  ir/l  with  both  Q(x,s)  and 
B(x,t),  and  none  of  Hx,s),  9(x,/)  and  9(x,c)  can  leave  t/  at  x.  A  similar  argument 
shows  that  refined  partition  edges  are  directed  away  from  dU  it  hitpoints. 

Corollary  2.9.3:  The  unrefined  Voronoi  diagram  V  forms  a  (directed)  tree  with  root 
c  and  edges  directed  towards  c. 

Proof:  No  cycles  are  possible  because  otherwise  some  Voronoi  cell  would  be 
separated  from  dU.    Only  c  has  out-degree  zero;  every  other  vertex  or  hitpoint  of  V 
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has  out-degree  1,  so  V  is  a  root-directed  tree.  D 

The  refined  Voronoi  diagram  V*  consists  of  V  together  with  refined  partition 
edges.  Each  such  edge  must  lie  entirely  within  a  single  Voronoi  cell,  having  one  end- 
point  on  dU  and  the  other  endpoint  at  a  reflex  comer  or  on  V  (by  Lemma  2.8.1). 

3.   The  Algorithm 

Procedure  gfv 

Input.  A  polygon  U  with  n  sides  and  a  set  SCU  of  k  sites. 

Output.  The  refined  furthest-site  geodesic  Voronoi  diagram  V* . 

1.  Triangulate  U. 

2.  Compute  C,  the  relative  convex  hull  of  S,  and  discard  all  non-extreme  sites  of  S. 

3.  Determine  two  or  three  two-fragment  instances  so  that  the  union  of  the  source  fragments  is  dl 

4.  Compute  V'CidUhy  calling  rgfs(u,v,s,t)  for  each  two-fragment  instance  iu,v,s,t). 

5.  Call  sweep  to  extend  V'  to  the  interior  of  U. 

Figure  3.1.1 


This  sertion  contains  the  algorithm  for  computing  V  .  An  outline  of  the  algo- 
rithm is  given  in  Figure  3.1.1.  Step  1  can  be  performed  in  time  0(/ilog/i)  [GJPT78]. 
The  relative  convex  hull  computation  of  the  second  step  can  be  accomplished  in  time 
0((n  +k) login  +k))  [T86].  The  third  step  also  takes  time  0((n  +k)  login  +k))  and  is 
described  in  Lemma  3.1.1  below.  The  fourth  step,  the  most  difficult  of  the  algo- 
rithm, is  the  computation  of  V*  restricted  to  dU.  It  is  discussed  in  Sections  3.1 
through  3.4.  The  last  step  is  the  extension  of  V*  to  the  interior  of  U.  This  is  done 
using  a  "reverse  geodesic  sweeping"  algorithm,  discussed  in  Section  3.5  below.  Both 
the  fourth  and  fifth  steps  take  time  Oi(n  +k) login  +k)). 

The  computation  of  V*  restriaed  to  dU  is  quite  similar  in  outline  to  Suri's  algo- 
rithm for  furthest  geodesic  neighbors  inside  a  simple  polygon[S87].  We  first  reduce 
the  computation  of  V*  on  dU  to  at  most  three  instances  of  the  "two-fragment  prob- 
lem." Roughly,  an  instance  of  the  two-fragment  problem  consists  of  a  fragment  of 
dU  a.nd  a  fragment  of  dC  so  that  all  furthest  sites  of  points  in  the  fragment  of  dU  are 
contained  in  the  fragment  of  dC.  Such  a  pair  must  also  satisfy  a  technical  condition 
given  below;  this  reduction  appears  in  Section  3.1.  The  algorithm  to  solve  the  two- 
fragment  problem  is  based  on  a  divide-and-conquer  schema  that  splits  an  instance 
into  two  smaller  instances.  The  basic  properties  of  the  divide-and-conquer  schema 
appear  in  Section  3.2.  Section  3.3  contains  the  exact  splitting  method  and  the  pro- 
cedures for  handling  the  base  cases  of  the  recursion.  The  complexity  analysis  appears 
in  Section  3.4.  We  show  that  the  sum  of  all  instance  sizes  at  each  level  of  recursion 
is  linear  in  n  +k.  This  implies  over  all  0((/i  +k)  login  +k))  running  time. 
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We  work  with  polygonal  relatively  convex  sets  in  addition  to  simple  polygons. 
Any  such  relatively  convex  set  Q  can  be  decomposed  into  a  collection  of  plateaus  and 
bridges.  Clearly  a  triangulation  of  Q  can  be  obtained  just  by  triangulating  each  pla- 
teau in  the  decomposition.  This  is  easily  done  in  time  0(m  logm)  [GJPT78],  if  m  is 
the  number  of  segments  in  the  boundary  of  Q.  Similarly  a  shortest  path  partition  of 
Q  from  an  arbitrary  point  in  it  can  be  obtained  by  using  a  shortest  path  partition 
algorithm  in  each  plateau  of  the  decomposition.  If  Q  has  been  triangulated,  this 
takes  time  0(m)  [GHLST]. 

3.1.  The  Two-Fragment  Problem 

A  two-fragment  instance  is  a  quadruple  (u,v,s,t)  where  u,v^dU,  m#v,  j  is  a  site 
furthest  from  u,  t  is  a.  site  furthest  from  v  (possibly  s=t),  and  g(u,v)  separates 
dU[u,v]  from  dC[s,t].  The  two-fragment  problem  is  "Given  two-fragment  instance 
(u,v,s,t),  compute  V*ndU[u,v]."  Observe  that  by  the  Ordering  Lemma,  only  the 
Voronoi  cells  of  sites  in  dC[s,t]  can  intersect  dU[u,v].  The  source  fragment  of  two- 
fragment  instance  (u,v,s,t)  is  3C/[m,v];  the  target  fragment  is  dC[s,t]. 

Lemma  3.1.1:  There  exists  a  set  of  at  most  three  instances  of  the  two-fragment  prob- 
lem so  that  the  union  of  the  source  fragments  is  dU.  The  instances  each  have  size 
0(n+k)  and  can  be  computed  in  time  0((n  4-A:)  log(n +/t))  given  a  triangulation  of 
U. 

Proof:  Choose  x,  y,  z,  w,  x* ,  y* ,  z*  as  in  the  proof  of  Lemma  2.7.6.  It  is  clear  that 
{x* ,y* ,y,z),  {y* ,z* ,z,w),  and  {z* ,x* ,w,y)  are  two-fragment  instances  each  of  size  at 
most  0{n-¥k)\  their  source  fragments  cover  dU.  As  for  computing  them,  the  choice 
of  X  was  arbitrary.  Site  y  can  be  determined  in  time  0((/i  +it)  log(/i  +/:))  by  comput- 
ing the  shortest  path  tree  from  x,  then  determining  the  distance  from  every  site  to  x 
using  a  planar  point  location  algorithm  in  the  resulting  shortest  path  partition.  Sites  z 
and  w  can  be  determined  similarly.  Surely  x* ,y* ,z*  can  be  computed  in  additional 
time  0(n).   The  case  when  j:=z  is  handled  similarly.   Q 

Let  D  be  the  relative  convex  hull  of  the  sites  on  flC[j,r].  Clearly  the  ordering 
of  sites  on  D  is  the  same  as  on  C,  with  the  addition  that  s  immediately  follows  t  in 
counterclockwise  order.  Let  R{,u.v,s,t)  be  the  relative  convex  hull  of  dU[u,v]  and  D. 
See  Figure  3.1.2. 

We  say  iu,v,s,t)  is  degenerate  if  D  is  contained  in  g{u,v).  If  {u,v,s,t)  is  degen- 
erate, there  can  be  no  sites  in  D  besides  s  and  t  and,  by  Lemma  2.7.2,  the  order 
along  g{u,v)  must  be  u,  t,  s,  v.   In  this  case  we  define  l  =  t  and  r=s. 

Suppose  (,u,v,s,t)  is  not  degenerate.  Let  /  be  the  counterclockwise  extreme  point 
of  D  from  u  and  r  the  clockwise  extreme  point  of  D  from  v.  By  Lemma  2.4.3, 
dU[u,v],  giv,r),  dD[r,l],  and  g{l,u)  constitute  a  counterclockwise  traversal  of  the 
boundary  of  R{u,v,s,t).  It  is  an  immediate  consequence  of  the  following  Lemma 
that  dD[r,l]  is  a  subpath  of  dC[s,t]. 

Lemma  3.1.2:  Sites  s,  r,  I,  t  are  in  that  counterclockwise  order  on  D,  not  necessarily 

all  distinct. 

Proof:  The  Lemma  is  trivial  if  r  =  l,  s  =  t,  or  if  {u,v,s,t)  is  degenerate,  so  assume  si^t 
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and  (u,v,s,t)  is  not  degenerate.  We  show  that  if  s  appears  on  dD  [r,l]  then  r=j  and 
if  r  appears  on  dD[r,l],  then  t=l.  Since  s  is  the  extreme  point  of  D  immediately 
counterclockwise  from  t,  this  implies  the  Lemma. 

Suppose  t  appears  on  dD[r,l]  and  f^l.  Since  s  is  immediately  counterclockwise 
from  /,  s  is  also  an  extreme  point  of  D  lying  on  dD  [r,l]  and  thus  an  extreme  point  of 
R(u,v,s,t).  Then  t,  s,  u,  and  v  appear  in  that  order  on  a  counterclockwise  traversal 
of  dR(u,v,s,t),  which  by  Lemma  2.3.7  implies  that  g(t,u)  meets  g(v,s),  contrary  to 
Lemma  2.7.2.  □ 

We  wish  to  give  a  definition  of  "left  side  connector"  and  "right  side  connector" 
to  capture  the  bounding  edges  of  R(u,v,s,t)  not  in  dC[r,l]  and  dU[u,v].  The  obvi- 
ous definitions  are  g{u,l)  and  g{v,r),  respectively.  Unfortunately,  these  definitions 
are  not  adequate.  In  Section  3.4,  we  analyze  the  size  of  side  connectors;  one  crucial 
property  used  in  our  argument  is  that  \i  s^t,  then  the  left  and  right  side  connectors 
are  disjoint  except  possibly  at  their  endpoints  (Lemma  3.1.3  below).  Unfortunately, 
this  is  not  true  for  side  connectors  defined  as  giu,l)  and  g(v,r).   See  Figure  3.1.2(b). 


(a)  (b) 

Figure  3.1.2.  Solid  outlines  D;  dashes  outline  R(u,v,s,t) 

It  is  clear  that  geodesic  giv,r)  has  connected  intersection  with  3D;  furthermore 
if  the  intersection  is  more  than  a  point  it  must  be  some  final  portion  of  geodesic 
gir'  ,r),  r'  the  site  of  D  immediately  clockwise  of  r.  Let  r  be  r  if  r=5,  otherwise  let  r 
be  the  first  point  of  g(r',r)  intersected  by  g(v,r).  The  right  side  connector  of 
(u,v,s,t)  is  g(v,r).  Similarly,  we  define  the  left  side  connector  of  iu,v,s,t)  to  be 
g{u,l)  where  /  is  /  if  l=t,  otherwise  /  is  the  first  point  of  gil',1)  intersected  by 
g{.u,l),  where  /'  is  the  site  of  D  immediately  counterclockwise  from  /.  See  Figure 
3.1.2.   A  connector  edge  is  a  link  in  either  the  left  or  right  connector. 

The  size  of  (u,v,j,r),  denoted  \{u,v,s,t)\,  is  |dC[5,r]|  + |3C/[u,v]|  plus  the  sizes 
of  the  side  connectors.  Since  dD[r,l'\  is  a  subpath  of  dC[s,t],  it  is  clear  that 
\dR(,u,v,s,t)\^\{u,v,s,t)\. 

Lemma  3.1.3:  If  the  side  connectors  of  a  two-fragment  instance  {u,v,s,t)  meet  at  a 
point  other  than  one  of  their  endpoints,  then  s  =  t. 
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Proof:  If  (u,v,s,t)  is  degenerate,  then  u,  l=t,  r=s,  v  appear  in  that  order  along 
giu,v),  and  the  side  connectors  are  disjoint  unless  /=r,^which  in  turn  forces  s  =  t.  So 
suppose  (u,v,s,t)  is  not  degenerate  and  point  a^u,v,r,l  is  common  to  both  side  con- 
nectors. Without  loss  of  generality  we  may  assume  that  a  is  a  reflex  comer  of  dU. 
Deleting  a  splits  R(,u,v,s,t)  into  two  components,  whose  closures  R\  and  Rizxe,  rela- 
tively convex  polygonal  regions  with  the  property  that  a  geodesic  connecting  two 
points  of  R{u,v,s,t)-{a)  passes  through  a  if  and  only  if  one  of  the  points  lies  in 
R\-{a}  and  the  other  in  R2-{a}.  Since  a  lies  on  the  left  connector  and  thus  on 
g(u,l),  either  u€/?i  and  /€/?2  or  vice  versa  (and  similarly  for  v  and  r).  Since 
dU[u,v],  g{v,r),  dC[r,l],  gil,u)  constitute  a  traversal  of  dR{u,v,s,t)  and  such  a 
traversal  cannot  visit  any  point  more  than  twice,  dU[u,v]  and  dC[r,l]  do  not  meet  a. 
In  particular,  u  and  v  lie  in  the  same  component.  Similarly,  r  and  /  must  lie  in  the 
same  component. 

Without  loss  of  generality,  assume  that  Ri  contains  u  and  v  and  Rj  contains  r 
and  /.  Now  also  siRi,  else  a  would  appear  on  dD[s,r]  and  thus  could  not  lie  on  a 
side  connector.  Hence  s€R2,  similarly  f€/?2.  Thus  s  is  furthest  from  m,  t  is  furthest 
from  V,  but  both  g(u,t)  and  g(v,s)  pass  through  a,  implying  s=t  by  Lemma  2.7.2.  a 

3.2.   The  Recursion  Scheme 

Procedure  rgfsiu,v,s,t) 

if  dU[u,v]  or  dC[s,t]  is  a  base  case 

then  compute  dU[u,v]nV*  directly 

else  choose  M'€3f/[M,v] 
locate  a  site  /(w)  furthest  from  w 
call  rgfs(u,w,s,f(w))  and  rgfsiw,v,f(,w),t) 

end 

Figure  3.2.1:  Recursive  procedure  rgfs. 

Figure  3.2.1  contains  a  recursive  procedure  rgfs  for  the  two-fragment  problem. 
Section  3.3  discusses  the  base  cases  and  the  choice  of  the  splitting  point  w,  while  the 
complexity  analysis  is  contained  in  Section  3.4.  We  now  discuss  some  basic  data 
structures  needed  for  the  recursion. 

At  each  level  of  recursion,  we  need  to  have  available  the  boundary  of  R(u,v,s,t) 
and  a  triangulation  of  its  interior.  For  the  topmost  level,  the  boundary  of  R(u,v,s,t) 
can  be  constructed  using  the  relative  convex  hull  algorithm  of  Toussaint[T86];  this 
takes  time  0((n +it)  log(n +ifc)).  Then  it  can  be  triangulated  in  additional  time 
0((rt+Jfc)log(/i+i!:))  [GJPT78]. 

For  the  recursive  step,  we  need  to  compute  /(w)  and  the  boundaries  and  tri- 
angulations  of  Riu,w,s,f{w))  and  R(w,v,f(w),t)  in  total  time  0{\{u,v,s,t)\).  To 
compute  /(w)  it  suffices  to  know  the  geodesic  distance  d„(r)  for  every  site  r  in 
dC[s,t];  d„(r)  can  be  determined  in  constant  time  if  the  cell  of  the  shortest  path 
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partition  of  R(u,v,s,t)  from  w  containing  r  is  known.  The  cell  containing  r  for  all 
sites  r  in  dC[s,t]  can  be  determined  in  total  time  Oi\(u,v,s,t)\)  as  follows.  We 
assume  the  shortest  path  partition  of  R(u,v,s,t)  from  w  has  been  computed  and 
refined  to  a  triangulation  (this  takes  only  linear  additional  time);  hence  each  cell  is  a 
triangle.  First  locate  the  triangle  containing  s;  this  clearly  can  be  done  in  the  allowed 
time  bound.  Then  traverse  dC[s,t],  in  one  step  moving  to  the  next  vertex  of  dC[s,t] 
or  to  the  next  intersection  of  the  current  edge  of  dC[s,t]  with  the  boundary  3A  of  the 
current  triangle  A  of  the  shortest  path  partition  from  w.  Notice  that  the  intersection 
of  dC  with  A  has  at  most  three  connected  components,  since  C  is  relatively  convex. 
Hence  the  traversal  of  dC[s,t]  takes  total  time  Oi\(u,v,s,t)\),  since  the  charge  for  a 
step  to  a  vertex  of  dC[s,t]  can  be  allotted  to  the  vertex  and  the  step  to  an  intersection 
with  A  can  be  allotted  to  one  of  the  three  connected  components  of  AfldC^.f]. 

We  compute  the  boundary  and  triangulation  of  R(u,w,s,f(w))  as  follows  (han- 
dling /?(w,v,/(w),f)  is  similar).  If  (u,w,s,fiw))  is  degenerate,  both  the  boundary 
and  triangulation  of  R(u,w,s,f(w))  can  be  easily  obtained  from  the  shortest  path  par- 
tition of  R(u,v,s,t)  from  w.  Otherwise  compute  r',  the  clockwise  extreme  point  of 
dC[s,f(w)]  from  w  using  the  shortest  path  partition  from  w.  Similarly  /',  the  coun- 
terclockwise extreme  point  of  dC[s,f{w)]  from  u,  can  be  determined  by  computing 
the  shortest  path  partition  of  R(u,v,s,t)  from  u.  Now  since  /'  is  extreme  in 
R(u,w,s,f(w)),  g(w,l')  splits  R(u,w,s,fiw))  into  two  pieces,  one  piece  lying  to  the 
left  and  one  to  the  right.  (Possibly  one  or  the  other  is  just  giw.l').)  A  triangulation 
of  the  piece  lying  to  the  right  can  be  obtained  by  refining  the  shortest  path  partition 
of  R(u,v,s,t)  from  w.  Similarly  a  triangulation  of  the  piece  lying  to  the  left  can  be 
obtained  by  refining  the  shortest  path  partition  of  R(u,v,s,t')  from  /'.  Notice  the 
links  in  giu,l'),  g(w,l')  and  g(w,r')  are  used  as  triangle  edges  in  this  triangulation. 
The  left  and  right  side  connectors  of  R(u,w,s,fiw))  are  easily  determined  from 
g(u,l')  and  giw,r').  This  computation  can  clearly  be  done  in  time  0(\iu,v,s,t)\). 

3.3.   Choosing  splitting  points  and  the  base  cases 

If  u  and  V  do  not  lie  on  the  same  wall  of  dU,  then  splitting  point  w  is  chosen 
simply  as  the  comer  of  dU[u,v]  so  that  |3C/[«,w]|  is  within  one  of  |aC/[H',v]|. 

If  dC[s,t]  is  a  single  site,  i.e.,  5=f,  then  all  of  dU[u,v]  lies  in  V(s).  We  need  to 
find  the  refined  partition  edges  of  V*  intersecting  dU[u,v];  it  is  sufficient  to  compute 
the  shortest  path  partition  of  Riu,v,s,t)  from  s,  which  can  be  done  in  time 
0{\(u,v,s,t)\)  given  the  triangulation  of /?(u,v,j,r)[GHLST]. 

It  is  possible  that  s^t  but  u  and  v  lie  on  the  same  wall  of  dU.  In  this  case  there 
is  no  obvious  splitting  point  w.  We  perform  a  "partition"  step:  segment  uv  is  split 
into  subsegments  so  that  within  each  subsegment  the  shortest  path  tree  from  any 
point  on  the  subsegment  to  the  sites  on  dC[s,t]  is  combinatorially  invariant.  This 
partitioning  is  described  below;  it  results  in  0{\(u,v,s,t)\)  subsegments  and  takes 
time  C>(  |(u,v,.y,r)  I  log|(M,v,.y,r)  I).  We  introduce  the  partition  points  as  dummy  ver- 
tices, and  use  them  as  splitting  points  in  the  divide  and  conquer.  Notice  that  the  par- 
titioning is  performed  at  most  once  on  a  path  from  the  topmost  instance  to  a  leaf 


29 


instance  in  the  recursion  tree.  Hence  we  introduce  only  0(n+k)  such  points,  since 
the  sum  of  instance  sizes  at  each  level  of  the  recursion  tree  is  0(n+k)  (see  the 
analysis  in  Section  3.4).  Similarly,  the  total  time  to  compute  the  partitioning  points  is 
0((n+k)\osin+k)). 

The  remaining  problem  is  to  handle  a  base  case  instance  (u,v,s,t),  where  s^^t 
and  the  shortest  path  tree  is  combinatorially  invariant  on  segment  uv=dU[u,v]. 
Because  of  this  invariance,  no  refined  partition  edge  of  V*  intersects  segment  uv,  in 
other  words  MvnV*=uvnV  is  the  set  of  bisector  hitpoints  on  uv.  Notice  that,  the 
partition  induced  on  segment  mv  by  V*  is  exactly  the  partition  induced  by  the  upper 
envelope  of  the  functions  dr,  where  r  is  a  site  in  dC[s,t].  Again  because  of  the  com- 
binatorial invariance  of  the  shortest  path  tree,  each  function  dr  is  "simple"  on  seg- 
ment Mv;  specifically  dr(x)  is  of  the  form  ci+vciix)  where  ci  is  constant  and  C2(x) 
depends  quadratically  upon  the  position  of  x  on  segment  mv.  (Observe  that  the  pur- 
pose of  partitioning  the  original  wall  was  to  ensure  that  ci  and  C2  are  fixed  over  the 
length  of  MV.  Their  values  for  each  site  r(:dC[s,t]  are  defined  by  the  identity  of  the 
anchor  of  x€mv  with  respect  to  r  and  the  distance  from  r  to  this  anchor;  all  of  this 
information  can  be  determined  in  linear  time  from  the  shortest-path  tree  of 
R(u,v,s,t)  from,  say,  m.)  Thus  in  constant  time  it  is  possible  to  determine,  for  a  pair 
of  sites  r  and  r'  the  (unique)  point  jc€mv,  if  any,  for  which  dr{x)=dr'(.x).  Now  by  the 
Ordering  Lemma,  Voronoi  cells  appear  along  segment  mv  in  the  same  order  as  the 
corresponding  sites  appear  along  dC[s,t].  This  implies  that  the  upper  envelope  of  the 
functions  dr  on  segment  mv  can  be  computed  in  time  proportional  to  the  number  of 
sites  (which  is  certainly  0(\(u,v,s,t)\)).  For  example,  an  incremental  algorithm  is 
sufficient.  Suppose  that  the  partition  of  segment  mv  induced  by  an  initial  subsequence 
of  the  sites  on  dC[s,t]  has  been  computed.  Then  the  partition  induced  by  adding  the 
next  site  in  order  can  be  determined  in  constant  time  plus  time  proportional  to  the 
number  of  cells  deleted  from  the  partition  of  segment  mv  computed  so  far. 

We  now  describe  the  partition  step,  which  uses  a  technique  similar  to  that  of 
[GHLST].  We  actually  partition  segment  mv  so  that  within  each  subsegment,  the 
shortest  path  tree  to  every  site  of  dC[s,t]  and  every  vertex  of  dRiu,v,s,t)  is  combina- 
torially invariant.  To  do  this,  compute  T(u),  the  shortest  path  tree  of  R(u,v,s,t) 
from  M.  Include  in  T(u)  the  links  of  geodesies  from  u  to  sites  of  dC[s,t]  not  appear- 
ing on  3/?(u,v,j,/)  (there  can  be  at  most  one  such  link  per  site  not  already  in  T{u)). 
The  augmented  ^(m)  can  be  computed  in  total  time  0(\(u,v,s,t)\)  using  the  tech- 
nique described  for  computation  of /(w)  in  Section  3.2.  Similarly  compute  T{v),  the 
augmented  shortest  path  tree  from  v. 

For  each  z,  z  either  a  site  of  dC[s,t]  or  a  vertex  of  dR(u,v,s,t),  compare  Q(z,u) 
with  6(z,v).  This  is  possible  using  T(u)  and  ^(v).  If  9(z,m)  =  0(z,v)  do  nothing.  If  z 
lies  on  segment  mv  do  nothing.  Otherwise  geodesic  triangle  Azmv  must  actually  be  a 
simple  polygon  (since  mv  is  a  wall).  Pass  a  line  through  the  first  link  of  ^(z,m);  it 
must  hit  segment  mv  since  the  interior  angles  of  Azmv  are  reflex  except  at  z,u,v.  The 
intersection  point  is  a  partition  point.  Similarly  obtain  a  partition  point  from  the  first 
link  of  giz,v).    It  is  easy  to  see  that  the  partition  points  generated  in  this  fashion 
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form  the  required  partition  of  segment  uv.  Generating  the  partition  points  takes 
0{\{u,v,s,t)  I)  time  and  sorting  them  requires  0( |(u,v,j,r) |  log  \{u,v,s,t)  |)  time. 

3.4.  Complexity  analysis 

Lemma  3.4.1:  Let  («,-,v,-,ji,r,)  be  a  topmost  two-fragment  instance.  There  are  at 
most  0(n  +k)  distinct  connector  edges  among  all  subinstances  of  (w,-,v;,5,-,r,). 
Proof:  By  the  discussion  of  the  preceding  section,  there  are,  over  the  course  of  algo- 
rithm execution,  only  0{n+k)  subinstances  of  the  two-fragment  problem,  hence  as 
many  connectors.  Since  ,?(m,-,v,)  separates  dU[ui,Vi]  from  3C[j,-,r,],  it  separates  the 
endpoints  of  each  connector  as  well.  Hence  by  Lemma  2.7.5,  there  are  at  most 
Oin+k)  distinct  connector  edges,  a 

Lemma  3.4.2:  Let  (M,-,v,-,j,-,r,)  be  a  topmost  two-fragment  instance.  At  each  level  of 
recursion,  each  connector  edge  appears  in  only  a  constant  number  of  subinstances  of 

(Ui,Vi,Si,ti). 

Proof:  We  count  the  number  of  times  edge  ab  can  appear  as  a  left  connector  edge 
directed  from  a  to  b.  Notice  a  and  b  may  be  assumed  to  be  reflex  comers  of  dU,  for 
otherwise  ab  is  necessarily  a  first  or  last  link  on  a  side  connector,  and  hence  can 
appear  in  only  one  left  connector  (three  if  ab  is  the  last  link  and  the  target  fragment 
consists  of  the  single  site  b  in  two  of  the  instances).  Let  Qa  be  the  set  of  points  x  of 
U  for  which  g{x,b)  passes  through  a.  Similarly,  define  Qb  as  the  set  of  points  that 
can  reach  a  only  through  b.  Clearly  Qa  and  Qb  are  non-empty  and  disjoint.  Let  A  be 
BU[ui,Vi]r\Qa  and  B  be  dC[si,ti]r\Qt,.  It  can  be  checked  that  A  is  a  single  fragment 
of  dU  and  that  B  is  a  single  fragment  of  dC.  Then  for  u€3t/[M,-,v,]  and  r€3C[5,-,r,], 
link  ab  appears  in  g(u,t)  exactly  if  m€A  and  t^B. 

We  first  claim  that  at  each  level  of  recursion  there  are  at  most  two  nonleaf 
instances  of  the  two-fragment  problem  for  which  both  the  source  fragment  interseas 
A  and  the  target  fragment  intersects  B.  To  see  this  suppose  (uj,Vj,Sj,tj),  j  =  1,2,3, 
are  instances  at  the  same  level,  dU[uj,Vj]nA=^0,  dC[sj,tj]r\B=^0,  and  ui,M2,M3  (and 
hence  si,S2,S3)  are  in  that  counterclockwise  order.  Since  these  instances  are  all  at 
the  same  level  of  recursion,  V2  appears  between  U2  and  M3  (possibly  V2=W3)  and  ^2 
appears  between  /i  and  t2  (possibly  ^2=^1)-  But  then  V2€A,  si^B,  so  the  geodesic 
g(v2,S2)  contains  link  ab.  Recalling  that  the  relative  interior  of  ab  does  not  lie  in  Qb 
and  hence  is  disjoint  from  dC[s2,t2]QB,  we  deduce  that  the  right  connector  of 
iu2,V2,S2,t2)  contains  the  relative  interior  of  ab.  Since  g(u2,t2)  also  contains  ab,  the 
relative  interior  of  ab  lies  in  the  left  connector  of  this  instance  as  well.  Hence  52=^2 
by  Lemma  3.1.3  and  (U2,V2,.S2.'2)  is  a  leaf  instance. 

We  now  claim  that  at  each  level  of  recursion  there  are  at  most  four  instances  of 
the  two-fragment  problem  (leaf  and  nonleaf)  with  ab  appearing  as  a  left  connector 
edge.  Such  an  instance  (u,v,s,t)  must  have  u  €A  and  r€5.  This  is  only  possible  if, 
for  its  parent  instance,  the  source  fragment  intersects  A  and  the  target  fragment  inter- 
sects B.   As  just  argued  there  are  only  two  such  parent  instances,  a 

Theorem  3.4.3:  V*ndU can  be  computed  in  time  0((n  +k)\og(n+k)). 

Proof:  Clearly  C/ can  be  triangulated  in  time  O(nlogn)  [GJPT78].    By  Lemma  3.1.1, 
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there  are  three  two-fragment  instances  with  union  of  source  fragments  equal  to  dU 
that  can  be  computed  in  time  0(n+k).  We  show  that  rgfs  solves  the  two-fragment 
problem  in  total  time  0((n  +k)  log(/i  +k))  for  each  (top-level)  instance,  proving  the 
Theorem. 

Consider  the  work  performed  by  rgfs  for  all  instances  at  a  particular  level  of 
recursion,  ignoring  recursive  calls  and  the  time  required  to  partition  walls  as  dis- 
cussed in  Section  3.3.  It  is  linear  in  instance  size  which  is  the  sum  of  the  sizes  of 
source  and  target  fragments  and  the  sizes  of  the  side  connectors.  The  total  size  of 
source  and  target  fragments  at  a  particular  level  of  recursion  is  0(n+k),  because 
source  and  target  fragments  are  partitioned  disjointly  except  for  endpoints,  and  there 
are  only  Oin+k)  possible  endpoints.  By  Lemmas  3.4.1  and  3.4.2,  the  total  size  of 
all  connectors  at  a  particular  level  of  recursion  is  also  0(n  +k).  Hence  the  total  work 
at  a  particular  level  of  recursion,  summed  over  all  instances  at  the  level,  is  0(n  +k). 

The  total  depth  of  recursion  is  0(\og(n+k)):  at  each  step  except  for  partition 
steps,  the  size  of  a  source  fragment  is  split  in  half.  At  a  partition  step,  the  size  of  the 
source  fragment  increases  to  at  most  0(n+k),  and  a  partition  step  happens  at  most 
once  on  a  path  in  the  recursion  tree  from  topmost  instance  to  leaf  instance. 

The  total  work  required  for  partitioning  is  0((n  +k)  log(n  +k)).  Hence  the  total 
work  to  solve  a  single  top-level  two-fragment  instance  is  0{(n  +k)  \og(n  +k)).   Q 

3.5.  Computiiig  V*. 

V*  is  computed  by  the  procedure  sweep  (Figure  3.5.1),  which  is  a  "reverse  geo- 
desic sweeping  algorithm";  it  progresses  from  5f/ towards  c,  the  center  of  C. 

Theorem  3.5.1:  Procedure  sweep  computes  V* .  It  can  be  implemented  to  run  in  time 
0(in +k)  log  in +k))  and  space  0(n+ifc). 

Proof:  For  positive  real  z,  let  Ddz)  be  the  geodesic  disc  of  radius  z  centered  at  c,  i.e., 
the  set  of  all  points  of  U  at  geodesic  distance  at  most  z  from  c.  We  claim  that  the 
while  loop  maintains  the  invariant  that  L  is  exactly  the  refined  Voronoi  edges  inter- 
sected by  dDcir),  in  order  around  dDdr).  This  follows  from  Lemma  2.9.2,  using 
standard  sweepline  arguments  [B079].  Hence  procedure  sweep  computes  V  . 

List  L  can  be  implemented  simply  as  a  circular  doubly-linked  list,  so  each  list 
operation  takes  constant  time.  Q  can  be  implemented  as  a  heap,  so  that  each  opera- 
tion takes  time  0(log(/i  +k)).  The  geodesic  center  c  of  C  can  be  computed  in  time 
0{{n +k)  log  (n+k))  as  follows.  Pollack,  Sharir,  and  Rote  [PSR87]  show  how  to 
compute  the  center  of  (the  set  of  vertices  of)  a  simple  polygon;  their  algorithm 
extends  easily  to  a  polygonal  relatively  convex  set  as  well.  Since  geodesies  restricted 
to  be  inside  a  relatively  convex  set  (with  respect  to  U)  are  identical  to  geodesies 
inside  U,  it  suffices  to  compute  the  center  of  C  =R(S). 

The  shortest  path  partition  of  U  from  c  can  be  computed  in  time  0(n)  since  U  is 
triangulated  [GHLST].  Given  the  shortest  path  partition  from  c,  the  geodesic  dis- 
tance from  a  point  x^U  to  c  can  be  computed  in  time  C>(logn),  using  a  planar  subdi- 
vision search  algorithm  to  locate  the  shortest  path  partition  cell  containing  x  (such  as 
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Input:  Triangulated  polygon  U,  refined  Voronoi  diagram  V*  restricted  to  dU. 

Output:  V' . 

Data  structures: 

L:  a  doubly-linked  circular  list  of  refined  Voronoi  edges. 

Q:  a  priority  queue  of  points  of  U,  ordered  by  decreasing  geodesic  distance  from  c. 

procedure  sweep 

Compute  c,  the  center  of  C. 
Compute  the  shortest  path  partition  of  U  from  c. 
Initialize  L  to  be  dUnv' . 

Initialize  Q  to  contain  all  pairwise  intersections  of  refined  Voronoi  edges  im- 
mediately adjacent  in  L  and  all  anchors  of  the  refined  partition  edges  appearing 
in  L. 
whaeG#0 

Extract  from  Q  the  point  v  of  maximum  geodesic  distance  from  c. 

Delete  from  L  all  refined  Voronoi  edges  with  head  v. 

Delete   from   Q  any   intersections  and/or  anchors  involving  just-deleted 

edges. 

if  V  is  not  an  anchor  then 
Insert  into  L  the  refined  bisector  edge  e  with  tail  v. 
Insert  into  Q  any  new  intersections  of  e  with  adjacent  refined  Voronoi 

edges. 

end  if 

end  while 
end  sweep 
Figure  3.5.1.   Procedure  sweep 

that  of  [ST]).  By  Theorem  3.4.3,  dUHV*  can  be  computed  in  time 
O  ((n +k)  login +k)).  Hence  L  and  Q  can  be  initialized  in  total  time 
Oan  +  k)  login +k)). 

Each  iteration  of  the  while  loop  uses  a  constant  number  of  operations  on  Q  and 
L  and  one  geodesic-distance  computation  for  each  item  inserted  in  Q;  hence  each 
iteration  takes  time  Oilogin+k)).  By  Lemma  2.8.3,  there  are  only  Oin+k)  itera- 
tions of  the  while  loop.  Thus  the  running  time  of  the  entire  algorithm  is 
Oiin+k)login+k)).  Clearly  the  space  usage  is  0(/i +/:).  □ 

Acknowledgement 

We  wish  to  thank  Micha  Sharir  for  suggesting  the  idea  of  the  proof  for  Lemma  2.5.3. 


-33 


References 

[A 87]  B.  Aronov,  "On  the  geodesic  Voronoi  diagram  of  point  sites  in  a  simple 
polygon,"  to  appear,  Algorithmica.  A  preliminary  version  appeared  in 
Proceedings  of  the  Third  Annual  Symposium  on  Computational  Geometry, 
1987,  pp.  39-49. 

[AT86]  T.  Asano,  G.T.  Toussaint,  Computing  the  Geodesic  Center  of  a  Simple 
Polygon,  in  Perspectives  in  Computing:  Discrete  Algorithms  and  Complexity, 
Proceedings  of  Japan-US  Joint  Seminar,  D.S.  Johnson,  A.  Nozaki,  T.  Nish- 
izeki,  H.  Willis,  eds,  June  1986,  pp.  65-79. 

[B079]  J.L.  Bentley,  T.A.  Ottmann,  "Algorithms  for  Reporting  and  Counting 
Geometric  Intersections",  IEEE  Transactions  on  Computers,  C-28,  pp.  643- 
647, (1979). 

[CG85]  B.  Chazelle,  L.  Guibas,  "Visibility  and  intersection  problems  in  plane 
geometry,"  Proceedings  of  the  Symposium  on  Computational  Geometry,  1985, 
pp.  135-146. 

[GHLST]L.  Guibas,  J.  Herschberger,  D.  Leven,  M.  Sharir,  R.E.  Tarjan,  "Linear 
time  algorithms  for  visibility  and  shortest  path  problems  inside  triangulated 
simple  polygons,"  Algorithmica  2,  1987,  pp.  209-233. 

[GJPT78]M.R.  Garey,  D.S.  Johnson,  F.P.  Preparata,  R.E.  Tarjan,  'Triangulating  a 
simple  polygon,"  Information  Processing  Letters,  Vol  7(4),  1978,  pp.  8-21. 

[LP84]  D.T.  Lee,  F.P.  Preparata,  "Euclidean  shortest  paths  in  the  presence  of  rec- 
tilinear barriers,"  Networks  14(3),  pp.  393-410. 

[PSR87]  R.  Pollack,  M.  Sharir,  G.  Rote,  "Computing  the  Geodesic  Center  of  a  Sim- 
ple Polygon,"  to  appear.  Discrete  and  Computational  Geometry. 

[PS85]  F.P.  Preparata,  M.L  Shamos,  Computational  Geometry:  an  Introduction, 
Springer-Verlag,  1985. 

[S87]  S.  Suri,  "Computing  geodesic  furthest  neighbors  in  simple  polygons,"  to 
appear.  Journal  Computer  and  Systems  Sciences.  A  preliminary  version 
appeared  in  Proceedings  of  the  Third  Annual  Symposium  on  Computational 
Geometry,  1987,  pp.  64-75. 

[ST]  N.  Samak,  R.E.  Tarjan,  'Tlanar  Point  Location  using  Persistent  Search 

Trees,"  CACM  Vol  29(7),  pp.  669-679,  July  1986. 

[T86]  G.  Toussaint,  "An  optimal  algorithm  for  computing  the  relative  convex  hull 
of  a  set  of  points  in  a  polygon,"  Signal  Processing  III:  Theories  and  Applica- 
tions. Proc.  qfEUSIPCO-86,  Part  2,  North-Holland,  1986,  pp.  853-856. 


NYU  COMPSCI  TR-420 
Aronov,  Boris 

The  furthest-site  geodesic 
Voronoi  diagram 

c.l 


NYU  COMPSCI  TR-420 
Aronov,  Boris 

The  furthest-site  geodesic 
Voronoi  diagram 

c.l 


DATE   DUE 


BORROWER'S  NAME 


This  book  may  be  kept 

FOURTEEN    DAYS 


A  fine  will  be  charged  for  each  day  the  book  is 

kept  overtinie. 

CAYLORD   142 

'H'NTED  IN  U    s    * 

